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A.  STATEMENT  OF  THE  PROBLEM  STUDIED 

The  work  carried  out  under  this  contract  consisted  of  theoretical  and 
computer  calculations  of  the  electrical  and  structural  properties  of  various 
point  defects  which  are  likely  to  be  electrically  active  in  CdTe  and  HgCdTe 
alloys.  These  electrically  active  defects  included  both  impurity  (indium, 
arsenic,  antimony)  and  self- interstitials  (  mercury,  cadmium,  tellurium).  The 
major  questions  addressed  in  this  work  were: 

1)  How  do  these  defects  diffuse,  and  what  determines  their  final 
distribution  in  the  lattice? 

2)  How  does  the  electrical  activity  of  these  defects  change  as  a  function 
of  their  position  in  the  lattice? 

Techniques  used  to  address  these  questions  included  a  tight-binding 

small-basis,  or  "local  matrix",  approach  for  estimating  point  defect  total 

energies,  electronic  levels,  and  the  character  of  the  localized  states, 

developed  for  covalently  bonded  materials  with  considerable  ionic  and  metallic 

character,  such  as  HgCdTe,'*'  ®  and  a  more  quantitatively  accurate  tight-binding 

supercell  approach,  which  was  used  to  investigate  the  effects  of  lattice 
7-8 

distortion.  The  local  matrix  method  was  primarily  developed  unde;  a 

previously  awarded  ARO  contract  (DAAG  29-85-0119),  and  it  is  more  fully 

discussed  in  the  reports  for  this  contract,  although  the  final  two  papers 

giving  an  extensive  description  of  this  method"*  ^  were  prepared  under  the 

present  contract.  The  bulk  of  the  work  supported  by  the  present  contract 

consisted  of  using  the  methods  developed  under  both  contracts  to  investigate 

problems  of  practical  interest  to  those  making  devices  using  CdTe  and  HgCdTe 

alloys,  such  as  the  electrical  inactivity  which  is  often  observed  for  a  large 

fraction  of  the  indium  added  to  CdTe  by  non-laser-assisted  molecular  beam 
9 

epitaxy  (MBE) . 
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This  contract  was  funded  in  part  by  the  Night  Vision  and  Electro-Optics 
Center  (NVEOC) ,  which  has  a  major  interest  in  developing  and  improving  devices 
based  on  CdTe  and  HgCdTe  alloys.  The  physical  insight  into  dopant  and  native 
defect  behavior  derived  from  the  theoretical  results  described  in  the  body  of 
this  report  was  made  available  to  researchers  working  with  Hg^  ^Cd^Te  and 
growing  InSb/CdTe  heteroj unctions  at  the  NVEOC  through  a  series  of  formal  and 
informal  talks. 

SIMS  experiments  on  diffusion  in  CdTe,  HgTe,  and  HgCdTe  alloys  at  the 
NVEOC,  designed  to  provide  data  for  theoretical  interpretation  and  to  test 
theoretical  conclusions,  were  planned  to  proceed  in  parallel  with  the 
theoretical  work  funded  under  this  contract.  Due  to  equipment  problems  and 
funding  constraints  at  NVEOC,  these  experiments  have  not  yet  progressed  beyond 
the  calibration  and  testing  stage,  and  comparison  of  the  theoretical  results 
obtained  with  experimental  results  for  diffusion  in  CdTe,  HgTe,  and  HgCdTe  was 

limited  to  a  review  of  the  existing  experimental  (and  theoretical) 

, .  10 
literature. 


B.  SUMMARY  OF  THE  MOST  IMPORTANT  RESULTS 


1.  Characterization  of  Interstitial  Defect  States 


11-13 

Previous  calculations  demonstrated  that  both  CdTe  and  Hg-rich  HgCdTe 

tend  to  have  high  total  concentrations  of  point  defects  during  growth  and 

annealing,  whether  or  not  electrically  active  impurities  have  been  added. 

Substantial  equilibrium  concentrations  of  both  vacancies  and  interstitials 

were  predicted  under  various  growth  and  annealing  conditions. 

11  13 

It  was  also  suggested  ’  that  the  macroscopic  quality  of  the  final 
crystal  depends  on  kinetic  processes  which  may  be  strongly  affected  by  lattice 
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strains  associated  with  impurities  such  as  zinc.  (These  strains  may  pin  the 

point  defects  generated  at  growth,  inhibiting  their  condensation  into  larger 

11  13 

extended  defects,  and  perhaps  explaining  *  the  dramatic  decrease  m 

14 

dislocation  densities  seen  in  CdQ  9gZnQ  Q^Te,  as  compared  with  CdTe.  ) 
Defect  diffusion  can  control  the  electrical  uniformity  and  stability  over  time 
of  doped  materials,  heteroj unctions ,  and  superlattices ,  as  well  as  the 
formation  of  macroscopic  defects  such  as  dislocations .  While  vacancies  and 
substitutional  and  antisite  defects  may  be  more  numerous,  interstitial  defects 
tend  to  move  easily  within  the  crystal  and  often  make  the  most  significant 
contribution  to  diffusion. 

Occupation  of  various  defect  sites,  ana  the  resulting  equilibrium  and 
kinetic  properties,  may  be  strongly  dependent  on  the  technological  processes 
used  to  prepare  the  alloy  for  device  applications.  For  example,  in 
nonequilibriura  situations  such  as  occur  in  material  processed  at  low 
temperatures,  quenched  after  high- temperature  growth  or  annealing,  or  damaged 
by  radiation,^  substitutional  impurities  may  be  preferentially  displaced  to 
interstitial  positions.  If  it  occurs,  this  displacement  will  affect  both  the 
diffusion  and  electrical  behavior  of  substitutional  dopants  such  as  indium. 

The  first  part  of  the  work  for  this  project  included  an  investigation  of 
the  total  energies  and  deep  electronic  defect  levels  associated  with  various 
impurity  and  self- interstitials  at  different  sites  in  the  lattice  (as  well  as 
a  few  substitutional  defects) ,  and  the  implications  of  these  results  for 
diffusion  and  electrical  behavior,  including  dopant  self- compensation.  The 
major  conclusions  are  summarized  in  Sections  a  and  b  below,  and  discussed  in 
greater  detail  in  the*  papers  included  as  Appendixes  A  and  B,  and  the  other 
cited  papers  included  in  the  list  of  publications. 

In  addition  to  accurately  showing  trends  in  the  electronic  defect  levels 
and  total  energies,  the  local  matrix  method  used  in  these  calculations  allows 
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us  to  form  a  good  qualitative  picture  of  the  charge  distribution  associated 
with  the  localized  defect  states.  In  the  absence  of  reproducible  experimental 
characterizations  of  the  localized  electronic  states  of  these  point  defects, 
the  calculated  wave  functions  may  be  used  to  estimate  carrier  capture  cross 
sections,  and  to  investigate  the  changes  in  bonding  at  the  defect  site  which 
may  occur  as  a  result  of  charge  carrier  capture.  The  local  matrix  method  also 
gives  a  particularly  clear  physical  understanding  of  the  role  of  the  local 
environment  in  determining  the  character  of  the  localized  defect  states.** 

As  an  example  of  this  characterization,  Figure  4  in  Appendix  D  shows  the 
formation  process  of  the  localized  defect  states  from  the  interstitial  atomic 
orbitals  and  the  sp -bonding  hybrids  of  the  neighboring  anions  and  cations  for 
the  mercury  interstitial  at  the  preferred  tetrahedral  positior  in  CdTe.  At 
each  stage,  the  fraction  of  the  total  charge  distribution  on  the  interstitial, 
the  nearest  neighbor  shell,  and  the  shell  of  atoms  bonded  to  the  nearest 
neighbors  is  given  for  each  localized  level.  The  wave  function  for  the  A^ 
tetrahedral  interstitial  level  in  the  gap  region  is  strongly  concentrated  on 
the  mercury  interstitial,  in  contrast  to  the  wave  function  for  the  A^ 
tetrahedral  interstitial  level  of  the  phosphorus  interstitial  (with  phosphorus 
nearest  neighbors)  in  GaP,  shown  in  Figure  1.  This  is  explained  as  a  result 
of  the  different  ordering  of  the  parent  levels  on  the  isolated  atoms:  both  the 
s  and  p  orbital  energies  for  the  mercury  interstitial  lie  above  the  energy  of 
the  sp-bonding  hybrids  on  the  nearest-neighbor  tellurium  atoms,  while  the  s 
orbital  energy  for  the  phosphorus  interstitial  lies  below  the  energy  of  the 
sp-bonding  hybrids  on  the  nearest-neighbor  phosphorus  atoms. 

a.  Electronic  Deep  Levels  and  Electrical  Activity 

In  order  to  test  the  reliability  of  the  local  matrix  method  for 
describing  the  localized  deep  and  hyperdeep  level  structure,  the  results  of 
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local  matrix  calculations  were  compared  with  self -consistent  Green's  function 

and  experimental  results  for  various  interstitials  in  Si  and  GaP,^  and  for 

2 

substitutional  defects,  including  vacancies,  in  Si,  GaP,  GaAs ,  and  CdTe.  The 

results  of  these  calculations  were  in  substantial  agreement  with  available 

experimental  data  for  cation  vacancies, ^  substitutional  gallium  and  indium  on 
17  17 

cation  sites,  and  bromine  on  the  anion  site  in  CdTe,  as  well  as  for  the 

various  defects  considered  in  Si,  GaP,  and  GaAs.  There  is  a  discrepancy  for 

oxygen  on  the  anion  site  in  CdTe,  since  a  level  in  the  gap  has  been  attributed 
18 

to  this  defect,  while  the  local  matrix  calculation  gives  an  level  for 

this  defect  higher  up,  at  2.2  eV  above  the  valence  band  edge,  in  the  absence 

of  lattice  distortion.  However,  this  discrepancy  may  be  due  to  the  existence 

of  substantial  lattice  distortion,  since  oxygen  is  a  small  impurity,  and  tends 

to  be  associated  with  lattice  distortion  in  other  systems. 

The  original  characterization  of  the  localized  defect  states  of  indium 

and  self- interstitials  (Hg,  Cd,  Te)  at  tetrahedral  and  hexagonal  positions  in 

H§1  xCVe  5iven  t^ie  PaPer  attached  as  Appendix  A.^  The  calculated 

energy  trends  and  conclusions  were  unchanged  by  a  subsequent  minor 

modification  of  the  parameters  and  the  inclusion  of  a  self-consistent 

5  6 

treatment  of  the  charge  distribution..  1  Two  major  conclusions,  discussed 

more  fully  in  Appendix  A,  emerged  from  this  work: 

(i)  A  level  of  A^  symmetry  due  to  the  mercury  interstitial  at  the 
tetrahedral  site  between  tellurium  nearest  neighbors  (which  was  assumed  to  be 
the  preferred  tetrahedral  site  for  mercury  interstitials)  was  identified  as  a 
possible  candidate  for  providing  the  midgap  tunneling  levels  in  Hg-rich 
HgCdTe . 

(ii)  It  was  shown  that  self- compensation  by  indium  interstitials  in  the 
preferred  tetrahedral  site  (assumed  to  be  between  tellurium  nearest  neighbors) 
cannot  account  for  the  electrical  inactivity  often  observed  for  a  large 
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fraction  of  the  indium  added  to  CdTe  by  non-photoassisted  MBE. 

Equilibrium  defect  densities  were  calculated  for  indium-doped  CdTe  at 

19 

temperatures  and  pressures  similar  to  those  which  occur  during  MBE  growth, 

11  13 

using  methods  previously  developed.  *  It  was  found  that  the  addition  of 

indium  promotes  the  creation  of  cation  vacancies.  At  low  temperatures  and 

high  doping  levels,  a  substantial  fraction  of  the  indium  is  bound  in  neutral 

complexes  consisting  of  a  single  cation  vacancy  and  two  substitutional  indium 

atoms,  while  nearly  equal  concentrations  of  positively  charged  substitutional 

indium  and  negatively  charged  complexes  of  a  cation  vacancy  with  a  single 

substitutional  indium  compensate  each  other,  leading  to  net  charge  carrier 

19 

concentrations  much  lower  than  the  total  indium  concentration. 

These  results  indicate  that  formation  of  indium-cation  vacancy  complexes 

can  fully  account  for  the  observed  electrical  inactivity,  if  the  material  is 

close  to  equilibrium  during  growth.  It  was  suggested  that  the  high  photon 

flux  occurring  during  photoassisted  MBE  could  excite  a  large  number  of  free 

carriers  at  the  growing  surface,  screening  the  Coulomb  attraction  between 

19 

charged  defects  and  inhibiting  the  formation  of  complexes.  X-ray  rocking 

studies  of  the  macroscopic  crystal  quality  of  samples  grown  by  photoassisted 

and  non-photoassisted  MBE  were  suggested,  to  help  determine  whether  the 

defects  are  able  to  approach  equilibrium  sufficiently  so  that  indium-cation 

vacancy  complex  formation  causes  self-compensation,  or  whether  the  inactive 

indium  is  bound  in  extended  (non-equilibrium)  defects  in  the  non-photoassisted 

samples.  In  the  latter  case,  it  was  suggested  that  laser  illumination  of  the 

substrate  or  the  impinging  atoms  could  provide  enough  energy  to  bring  the 

system  closer  to  equilibrium,  and  reduce  concentrations  of  these 

19 

non- equilibrium  defects. 
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b.  Total  Energies  and  Diffusion  Barriers 

The  local  matrix  results  for  the  total  energies  of  interstitial  defects 

in  various  positions  in  Hg-  Cd  Te  are  presented  in  the  paper  attached  as 
3 

Appendix  B.  The  major  conclusions  of  this  work,  described  in  more  detail  in 
Appendix  B,  are  summarized  below: 

(i)  The  preferred  tetrahedral  interstitial  site  for  mercury,  cadmium,  and 
indium  in  Hg1  Cd  Te  is  between  tellurium  nearest  neighbors.  These  results 
support  the  earlier  assumptions  we  made^  in  identifying  the  mercury 
interstitial  at  this  site  as  a  likely  candidate  for  providing  midgap  tunneling 
levels  in  Hg-rich  HgCdTe,  and  arguing  that  self-compensation  arising  from 
indi  an  interstitials  is  unlikely. 

(ii)  The  preferred  tetrahedral  interstitial  site  for  tellurium  in 
Hg1  Cd  Te  is  between  cation  nearest  neighbors. 

(iii)  Charge  carrier  capture  lowers  diffusion  barriers  for  simple 
interstitial  diffusion  down  the  open  channels  in  the  lattice  for  all  the 
interstitials  considered. 

(iv)  Diffusion  down  these  open  channels  may  be  important  for  mercury 
interstitials,  but  it  appears  to  be  less  favorable  for  indium,  and  very 
unfavorable  for  cadmium,  which  must  diffuse  by  some  more  complicated 
mechanism,  probably  involving  exchange  with  atoms  on  the  lattice  sites. 
(These  conclusions  were  in  agreement  with  simultaneously  reported  experimental 
results. 

2.  Effects  of  Lattice  Distortion  on  Defect  Properties  and  Interactions 

Self-consistent  tight-binding  supercell  calculations  were  carried  out  to 
investigate  the  lattice  relaxation  associated  with  vacancies  and  various 
impurity  (indium,  arsenic,  antimony)  and  self- interstitial  (mercury,  cadmium, 
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tellurium)  defects  in  Hg^  ^Cd^Te.  *  The  results  of  this  work  are  summarized 

g 

in  the  paper  attached  as  Appendix  C.  The  major  conclusions  drawn  from  this 
investigation,  discussed  in  more  detail  in  Appendix  C,  are  listed  below: 

(i)  The  strain-mediated  attraction  found  between  tellurium  xnterstitials 
and  cation  vacancies  may  promote  the  creation  of  antisite  defects. 

(ii)  The  large  strain  found  for  tellurium  interstitials  may  help  explain 
the  tendency  for  excess  tellurium  to  condense  into  tellurium  inclusions. 

(iii)  The  unfavorability  of  simple  interstitial  diffusion  of  cadmium  down 

3 

the  open  channels  in  the  lattice  which  was  predicted  theoretically  and 

20 

verified  experimentally  is  seen  to  be  a  consequence  of  the  large  strain 

energy  for  cadmium  interstitials  in  the  hexagonal  and  alternate  tetrahedral 
sites. 

(iv)  The  localized  A^  electronic  state  of  the  tetrahedral  mercury 

interstitial  identified  above^*  becomes  slightly  more  localized  on  the 

interstitial  as  a  result  of  relaxation  of  the  surrounding  lattice,  but  it 

remains  a  candidate  for  providing  midgap  tunneling  levels  in  Hg-rich  HgCdTe. 

Further  work  on  improving  the  accuracy  of  these  tight-binding  supercell 

total  energy  calculations  for  a  wide  range  of  defects  is  still  in  progress.  As 

a  part  of  this  work,  a  preliminary  investigation  of  the  role  of  the  d  states 

21 

in  determining  the  behavior  of  defects  in  CdTe  and  related  alloys  has  begun. 
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REPORT  OF  INVENTIONS 


Local  Matrix  Method: 

A  new  theoretical  approach  for  characterizing  the  localized  defect  states 
due  to  microscopic  structural  defects  was  developed,  tested,  and  used  to  study 
defects  in  Hg..  Cd  Te  under  the  present  contract  (DAAL  03-87-K-0061)  and  a 
previously  awarded  ARO  contract  (DAAG  29-85-K-0119) .  This  method  and  its 
applications  are  described  briefly  in  the  final  report  for  DAAG  29-85-K-0119 , 
in  the  papers  attached  to  this  report  as  Appendixes  A,  B,  and  D,  and  in  the 
other  publications  cited  as  references  in  this  report.  An  extensive 
presentation  of  the  method  is  given  in  Goettig  and  Morgan-Pond,  Phy.  Rev.  B 
42,  11730  (1990). 
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APPENDIX  A  A-i 

Deep  interstitial  levels  jn  Hg1_xCdxTe 

S.  Goettiga)  andC.G.  Morgan-Pond 

Department  of  Physics  and  Astronomy,  Wayne  State  University,  Detroit ,  Michigan  48202 
( Received  3  November  1987;  accepted  29  February  1988) 

The  deep  levels  of  indium  and  self-interstitials  ( Hg,Cd,Te)  in  tetrahedral  and  hexagonal  positions 
in  Hg,  _  xCdtTe  are  characterized  and  energy  trends  are  calculated  within  the  defect-molecule 
approach.  We  have  extended  this  method  to  include  metallic  and  polar,  as  well  as  covalent, 
coupling.  Possibilities  for  self-compensation  due  to  interstitials  in  In-doped  He,  _  xCdxTe,  and  a 
possible  candidate  for  the  midgap  tunneling  centers  in  Hgo.8Cd0  2Te  are  discussed. 


I.  INTRODUCTION 

The  identification  of  deep  levels  of  various  defects  in 
Hg,  _  .rCdxTe  is  an  important  step  towards  the  understand¬ 
ing  of  the  defect-induced  behavior  of  these  materials.  While 
most  of  the  unambiguously  identified  levels  have  been  attri¬ 
buted  to  substitutional  impurities,1  knowledge  about  inter¬ 
stitial  defects  in  HgCdTe  is  still  rather  fragmentary. 

In  this  paper  we  develop  a  defect-molecule  theory  of  deep 
localized  states  of  neutral  interstitial  defects  in 
Hg,_vCdtTe.  The  linear  combination  of  atomic  orbitals 
( LCAO )  defect-molecule  method  allows  us  to  calculate  the 
energies  of  the  deep  levels,  their  symmetries  and  occupan¬ 
cies,  as  well  as  the  corresponding  electronic  charge  densities. 
We  extend  the  method,5  recently  applied  to  self-interstitials 
in  Si,6  to  account  for  the  metallic  coupling  ( promotion  ener¬ 
gy)  and  the  polar  coupling,  both  of  which  an  important  in 
HgCdTe.  This  simple  model  allows  us  to  form  a  qualitative 
picture  of  the  localized  interstitial  levels,  and  to  investigate 
energy  trends  for  the  deep  levels.  Ongoing  work16  has  con¬ 
firmed  the  insensitivity  of  the  trends  noted  here  to  some  of 
the  approximations  made.  However,  the  calculated  energy 
levels  cannot  be  viewed  as  quantitatively  reliable. 

We  consider  a  molecule  (cluster)  composed  of  the  inter¬ 
stitial,  its  nearest  neighbors,  and  all  atoms  bonded  to  the 
nearest  neighbors.  The  cluster  is  separated  from  the  bulk  by 
systematically  neglecting  the  promotion  energy  on  all  atoms 
bonded  to  the  nearest  neighbors  of  the  interstitial,  with  no 
additional  boundary  condition  imposed.  Both  tetrahedral 
and  hexagonal  positions  of  the  interstitial  are  considered. 
The  position  of  the  valence  band  is  calculated  using  the  full 
sp 3  tight-binding  scheme  with  the  same  parameters  as  are 
used  in  the  defect  level  calculations. 

In  applying  the  formalism  to  HgCdTe,  we  have  concen¬ 
trated  on  indium  and  self-interstitials.  Self-interstitials  are  of 
general  interest  as  potentially  ever-present  elementary  de¬ 
fects.  It  has  also  been  suggested2  that  they  might  provide 
deep  levels  for  the  midgap  tunneling  observed  in 
HgasCdo^Te.  The  identification  of  such  levels  may  have 
important  practical  applications. 

The  doping  of  Hg^^Cd^Te  with  indium  has  attracted 
much  experimental  attention,  but  is  still  not  fully  under¬ 
stood.  When  introduced  in  small  concentrations  at  high  tem¬ 
perature,  indium  acts  primarily  as  a  substitutional  dopant 
with  shallow  donor  levels.  However,  when  indium  is  added 
to  CdTe  films  during  low-temperature  [molecular-beam 
epitaxy  (MBE)]  growth,  a  large  fraction  of  the  indium  is 


observed  to  be  electrically  inactive  (except  when  added  by 
photoassisted  MBE3).  It  is  known  that  substitutional  impur¬ 
ities  may  be  preferentially  displaced  to  interstitial  sites  in 
material  with  high  (nonequilibrium)  densities  of  native  de¬ 
fects  resulting  from  radiation  damage4  or  low-temperature 
growth.  Thus,  the  electrical  inactivity  of  indium  when  added 
to  CdTe  during  low-temperature  growth  may  be  attributed 
to  preferential  displacement  of  In  atoms  from  substitutional 
to  interstitial  sites,  resulting  in  formation  of  interstitial  deep 
acceptor  levels  which  lead  to  self-compensation.  We  discuss 
this  explanation  based  on  the  indium  interstitial  deep  levels 
calculated  within  the  present  model. 

II.  THEORY 

A.  Defect-molecule  model 

Let  us  consider  a  tetrahedrally  bonded  crystal  composed 
of /V  atoms  of  P  different  kinds  ( a=  1,2, ...,/>)  and  a  neutral 
interstitial  defect  (a  —  d).  We  assume  that  both  host  and 
interstitial  atoms  have  only  s  and  p  valence  electrons,  i.e.,  all 
d  shells  are  closed  and  do  not  take  part  in  the  bonding.  As  a 
basis  for  the  LCAO  expansion,  we  can  use  one  5-like  ($fl)) 
and  three  p-like  (<f>(xa\<f>ya\<f>{za))  quasiatomic  orbitals  cen¬ 
tered  on  each  host  and  defect  site.  Equivalently,  we  may  use 
sp*  tetrahedrally  directed  hybrids  <p/(/,  where  denotes  a 
hybrid  centered  on  the  /th  site  with  the  positive  lobe  pointing 
towards  the yth  nearest-neighbor  site.  Thus,  the  one-electron 
wave  function  of  the  system  can  be  written  in  the  convenient 
LCAO  form: 

*(«•)  =  2  I  <XU<Pl.J^  +  X  «n^rf>(r)  *  0) 

I  jU)  «■  1 

with  the  first  term  referring  to  the  host  and  the  second  to  the 
defect.  The  origin  of  the  reference  frame  was  chosen  at  the 
interstitial  site.  The  first  index  in  cpitJ  is  multiple  and  de¬ 
scribes  both  the  position  of  the  site  r,  and  the  kind  of  atom 
occupying  this  site.  The  summation  over  j  in  Eq.  ( 1 )  is  over 
the  four  nearest-neighbor  sites  of  the  /th  site.  We  use  the 
Schrodinger  equation  HV  =  EV,  where  H  is  the  total  one- 
electron  Hamiltonian  consisting  of  the  kinetic  part  and  the 
potential  part  F(r),  due  to  interaction  with  the  atomic  cores 
and  the  ( self-consistent )  field  due  to  all  other  electrons.  This 
gives  the  matrix  equation 

X  X  auj^<Pr./W\(Ptj)  ESlfSjf}  =  0f  (2) 

imOJU)  m  1 

where  the  compact  notation  tp0)rt  a0i„  has  been 
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introduced.  All  orbital  overlaps  were  neglected  in  Eq.  (2). 
Considering  on-site  matrix  elements  in  Eq.  (2),  for  th«- 


host  sites  (/#0)  we  have 

/  .m  v  f(l/4)(£‘a,  +  3£^)=£a 
(<pj  k./>  1  (i/4)(£‘«>  -E\a))  =  Aa 

for;-/ 

for#/  ’ 
(3) 

where 

£r  =  <rwr>. 

(4a) 

wr>* 

(4b) 

lated  energy  levels  often  turn  out  empirically9  to  be  closer  to 
the  true  macroscopic  limit  for  clusters  which  include  only 
bonds  involving  under-  or  overcoordinated  atoms  than  for 
clusters  which  include  a  moderate  number  of  bonds  between 
properly  coordinated  atoms  as  well.  Although  the  actual 
charge  density  associated  with  a  defect  level  extends  beyond 
the  boundaries  of  the  small  cluster  we  have  chosen,  the  rela¬ 
tive  charge  densities  on  the  different  ions  closest  to  the  defect 
have  been  shown  to  be  relatively  unaffected  in  calculations  of 
this  kind  when  the  size  of  the  cluster  is  extended.^ 


Polar  coupling  (in  nonhomopolar  crystals)  is  due  to  differ¬ 
ent  values  of  Ea  for  anions  and  cations,  while  metallic  cou¬ 
pling  is  due  to  nonzero  promotion  energy  4AU .  We  are  as¬ 
suming  here  that  the  on-site  matrix  elements  are  determined 
by  the  local  (on-site)  behavior  of  the  potential  F(r),  neglect¬ 
ing  the  influence  of  the  interstitial  on  their  values.  ( Ea  and 
Ad  depend  on  the  kind  of  atom  but  not  its  position  relative  to 
the  interstitial.)  The  nonvanishing  on-site  matrix  elements 
for  the  defect  site  are  given  by  Eq.  (4)  with  a  =  d. 

The  off-site  matrix  elements  in  Eq.  (2)  are  considered 
within  the  nearest-neighbor  and  two-center7  approxima¬ 
tions,  i.e.,  they  are  described  by  four  parameters:  VSi(J ,  Vspa , 
Vppai  and  Vpp7T.  Considering  matrix  elements  between  host 
hybrids,  we  make  the  further  simplification  of  including  only 
the  interaction  between  hybrids  pointing  towards  each  other 
(covalent  coupling) 

(<PiJH\<pJtl)  =  /?,  (5) 

neglecting  here  other,  weaker,5  interactions  between  near¬ 
est-neighbor  hybrids.  Here  fi  =  ( V^a  -  2V$Vspo 
—  3Kp/w)/4.  The  interaction  between  the  defect  and  its 
nearest  neighbors  will  be  considered  in  the  next  subsections. 

The  separation  of  the  defect  molecule  from  the  crystal  is 
an  approximate  procedure  in  which  certain  matrix  elements 
in  Eq.  (2)  are  neglected.  We  have  chosen  the  molecule  to 
consist  of  the  interstitital,  its  first  nearest  neighbors  ( 1  nn)t 
and  all  atoms  bonded  to  first  nearest  neighbors  ( 2nn ),  in 
both  tetrahedral  and  hexagonal  positions  (Fig.  1).  This  is 
the  smallest  reasonable  cluster  which  can  be  chosen  without 
breaking  any  covalent  bonds  (the  strongest  coupling  in  the 
crystal).  The  separation  of  this  molecule  from  the  crystal 
consists  of  neglecting  Afl  (the  promotion  energy)  on  all  Inn 
atoms. 

The  choice  of  the  cluster  is  important,  and  adding  a  few 
more  shells  to  a  small  cluster  does  not  always  give  better 
results.8  In  choosing  the  cluster,  we  were  led  by  two  criteria. 
First,  all  bonds  which  are  most  affected  by  the  presence  of 
the  interstitial  (i.e.,  bonding  involving  1  nn  atoms,  which  are 
overcoordinated  due  to  the  presence  of  the  interstitial)  were 
systematically  included.  Second,  inclusion  of  bulklike  bonds 
between  tetrahedrally  coordinated  host  atoms  was  avoided. 
These  bonds  are  included  in  calculations  for  large  clusters, 
consisting  of 2000  or  more  atoms.  However,  inclusion  of  less 
than  a  macroscopically  large  number  of  bulklike  bonds  in 
small-cluster  calculations  causes  the  mixing  in  of  poorly  de¬ 
scribed  extended  states  which  are  not  true  bulk  states.  Al¬ 
though  this  correctly  represents  the  tendency  of  the  charge 
distribution  to  spread  out  beyond  the  Inn  atoms,  the  calcu- 


B.  Tetrahedral  interstitial 

In  this  case,  the  molecule  consists  of  15  atoms  and  in¬ 
volves  36  basis  states.  Four  states  are  located  on  the  intersti¬ 
tial.  There  are  four  Inn  atoms,  each  with  a  total  of  eight 
states  involved  in  its  four  bonds  to  Inn  atoms,  thus  account¬ 
ing  for  the  32  other  basis  states.  Since  six  of  the  Inn  atoms  are 
bonded  to  two  1  nn  atoms,  there  are  ten  Inn  atoms  in  all  [see 
Fig.  1(a)).  We  note  that  because  Atf  on  Inn  atoms  was  ne¬ 
glected,  the  molecule  can  be  seen  as  consisting  of  the  intersti¬ 
tial  and  four  side  clusters  (each  containing  five  host  atoms) 
coupled  to  each  other  only  through  the  interstitial. 


Fig.  1.  Defect  molecules  used  for  (a)  tetrahedral  and  (b)  hexagonal  inter¬ 
stitials  in  the  zinc-blende  lattice.  The  defect  is  indicated  by  «.  the  I  nn 
atoms  are  marked  with  a  horizontal  bar,  and  2 nn  atoms  are  marked  with  a 
vertical  bar.  Cations  and  anions  are  indicated  by  O  and  •  ( or  vice  versa ) 
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Since  the  interstitial  is  tetrahedrally  coordinated,  it  is  con¬ 
venient  to  choose  the  basis  on  the  defect  to  be  four  tetrahe¬ 
drally  directed  hybrids.  Then  we  can  assume  that  the  hy¬ 
brids  on  a  given  1/m  atom  couple  only  to  the  interstitial 
hybrid  which  points  toward  this  atom.  The  couplings  with 
the  other  three  interstitial  hybrids  are  roughly  an  order  of 
magnitude  smaller,  and  shall  not  be  considered  here.  The 
addition  of  these  couplings  changes  the  energies  of  the  deep 
levels  within  1  eV  of  the  gap  by  <0.2  eV,l(> 

The  36x36  matrix  for  the  tetrahedral  interstitial  in  the 
zinc-blende  crystal  (P=  2)  was  diagonalized  analytically 
Omitting  the  details  of  the  calculation,  we  only  mention  that 
the  diagonalization  relies  on  representing  the  Hamiltonian 
matrix  in  Eq.  ( 2 )  in  the  basis  consisting  of  the  eigenfunctions 
of  the  five-atom  side  clusters,  and  then  using  the  Lowdin10 
transformation.  The  deep-level  energies  £  are  then  given  ex¬ 
actly  as  the  solutions  of  the  following  fifth-order  equations: 

£-tfJJSdj£i  +  jSiJ*l 

£-£s,  £-£,_ 


4-  3 


!cjj- 


"Pi  ,  K- 


\*P 


(6) 


E-Ep+  E  —  Ep. 

with  €  =  E[d>  for  the  nondegenerate  state  of  <4,  symmetry, 
ande  =  Epd>  for  the  threefold  degenerate  state  of  ^symme¬ 
try.  Here, 

£J±  =  ( 1/2)  {£,  —  3  A ,  +  £, 


±V(£,-3A  ,-£,); +  4/^},  (7) 

Ep±  =  ( 1/2)  {£,  +  A,  +  £, 


±  n/(£,  +  A,  -£,)-  +  W),  (8) 

cJ±  =(l/2){l+£V(£ti  -£2)2}-|/2,  (9) 

cp±  *(1/2){1  +02/{Ep±  —  £2)2}~i/:,  (10) 

are  the  eigenvalues  and  the  normalization  factors  for  the 
eigenfunctions  of  the  five-atom  cluster.  We  note  that  all  1/m 
atoms  are  of  the  same  kind  ( either  cations  or  anions ) ,  and  all 
2/m  atoms  are  of  the  opposite  kind.  In  Eqs.  (7)  and  (8),  ft 
and  A,  refer  to  1/m  atoms,  while  E2  refers  to  2nn  atoms.  The 
parameters  ft  =  Vxiif  -  3  ft,,,  and  ft  =  -  Kpa  -  Vppo 
describe  the  coupling  of  the  s  and  p  orbitals  on  the  Inn  atom. 
We  note  that  in  the  absence  of  the  interstitial  (e-  oo ),  the 
energy  levels  of  the  defect  molecule  are  given  by  Eqs.  ( 7 )  and 
(8). 

The  formulas  obtained  for  the  corresponding  wave  func¬ 
tions  are  rather  complicated  and  will  not  be  presented  here. 


C.  Hexagonal  interstitial 

In  this  case  the  molecule  consists  of  19  atoms  and  involves 
40  basis  states.  Again,  four  of  these  states  are  located  on  the 
interstitial.  Each  of  the  six  1/m  atoms  is  bonded  to  two  2 nn 
atoms  and  two  other  Inn  atoms,  so  that  the  six  Inn  atoms 
form  a  puckered  ring  around  the  interstitial  [see  Fig.  1(b)]. 
Each  unit  consisting  of  one  1  nn  atom  and  the  two  2/m  atoms 
it  bonds  to  contains  two  complete  bonds  and  the  two  spy 
orbitals  pointing  toward  the  neighboring  Inn  atoms.  Thus, 
we  need  a  total  of  six  basis  states  per  unit,  or  36  basis  states 
for  all  six  units. 
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The  interstitial  defect  now  has  lower  symmetry  than  in  the 
tetrahedral  case.  The  coupling  of  the  interstitial  with  its 
nearest  neighbors  was  considered  exactly,  and  expressed  in 
terms  of  all  four  interactions  ft  m.  The  40x40  matrix  was 
diagonalized  numerically.  The  resulting  levels  are  either  sin¬ 
gly  degenerate  (A)  or  doubly  degenerate  (£). 

D.  Choice  of  parameters 

The  on-site  energies  E  and  E pa)  were  determined  from 
atomic  levels.  The  s  energies  for  Hg  and  Cd  and  the  p  ener¬ 
gies  for  Te  and  In  were  obtained  from  experimental  ioniza¬ 
tion  energies, u  by  averaging  over  all  states  with  the  ground- 
state  configuration  (r,  s2pA ,  etc.).  The  other  energies  were 
obtained  from  the  above  energies  by  using  estimates  of  the 
s-p  promotion  energy  based  on  Hartree-Fock  calcula¬ 
tions.12  The  resulting  energies  were  ft  =  —  10.43,  —  8.99, 
-  18.0,  and  -  1 1.4  eV  and  ft  =  -  5.0,  -  4.77,  -  7.94, 
and  —  5.64  eV  for  Hg,  Cd,  Te,  and  In,  respectively. 

The  off-site  parameters  were  assumed  to  depend  only  on 
the  bond  length  ft  and  were  determined  by  VlVm  =  i)lVm  fr/ 
tnd 2  with  the  tight-binding  values  ?/„„  =  —  1.32, 
VsPa  =  1-42,  Vppo  =  2.22,  and  Vppn  =  -  0.63. n 

III.  RESULTS  AND  DISCUSSION 

We  obtained  the  energies  of  deep  levels  in  the  energy-gap 
region  due  to  In  and  self-interstitials  (Hg,  Cd,  Te)  for  CdTe 
and  HgTe.  These  are  shown  in  Fig.  2.  The  position  of  the 
defect  in  the  lattice  is  given  in  parentheses.  The  occupancy  of 
the  neutral  interstitial  levels,  their  symmetries,  and  the  cor¬ 
responding  one-electron  charge  densities  on  the  defect,  on  all 
1  ,  atoms,  and  all  2  nn  atoms  are  also  indicated. 

The  position  of  the  valence  band  (and  the  energy  gap 
shown  for  CdTe)  was  determined  by  the  full sp*  tight-bind¬ 
ing  theory,  using  the  same  values  of  VWm  as  in  the  interstitial 
calculation.  In  Fig.  2,  we  have  shown  the  defect  level  ener¬ 
gies  relative  to  the  calculated  valence-band  edge  for  CdTe 
and  HgTe.  Use  of  the  valence-band  edge  as  the  reference 
level,  rather  than  the  conduction-band  edge,  appears  to  give 
good  results  for  the  deep  levels  of  the  vacancy  in  several 
materials,  as  will  be  discussed  elsewhere. 16  Although  the  cal¬ 
culated  energy  gap  does  not  vanish  in  HgTe,  due  to  the  ne¬ 
glect  of  the  spin-orbit  interaction,  the  gap  in  Fig.  2  is  drawn 
to  show  the  energy  gap  going  to  zero  as  the  alloy  becomes 
more  Hg  rich,  as  is  seen  experimentally.  The  lines  between 
the  localized  defect  levels  show  schematically  the  behavior 
of  the  levels  as  we  go  from  CdTe  to  HgTe. 

The  levels  were  analyzed  with  respect  to  their  localized 
deep  character,  to  avoid  inclusion  of  quasiextended  (or  sur¬ 
face)  states.  In  this  analysis  the  following  criteria  were  used 
to  identify  localized  states.  ( i )  large  shift  of  the  energy  of  the 
level  as  compared  to  bulk,  or  isolated-defect  parent  level, 
(ii)  large  amount  of  charge  on  the  defect  site,  and  (iii)  sig¬ 
nificant  change  of  wave  function  on  \nn  and  2 nn  atoms  as 
compared  with  the  parent  bulk  state.  Using  these  criteria,  we 
find  that  the  hexagonal  level  of  Hg  at  £  =  ft  +1.2  eV  in 
HgTe  loses  its  localized  character  in  CdTe,  and  is  thus  not 
shown  on  the  CdTe  side.  Similarly,  the  hexagonal  level  of  Hg 
at  £  =  ft.  +  0.9  eV  in  CdTe  was  not  recognized  as  localized 
on  HgTe  side. 
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CdTe  - -  HgTe 


Fig  2.  The  localized  levels  due  to  In  and  self-interstitials  in  tetrahedral  and  hexagonal  positions  tn  HgCdTe.  X(b.  Yh  tetrahedral  level  of  X interstitial  located 
between  Y atoms;  (hex):  hexagonal  level;  A,  E,A,t  7V  symmetry  of  the  state  —OQOO  .  occupany  of  the  state;  :  occupied;  *0* ;  empty,  (a.&c): 
charge  densities  on  a,  interstitial;  b,  all  1  nn  atoms;  c,  all  Inti  atoms. 


Considering  the  behavior  of  the  deep  levels  and  their  posi¬ 
tions  in  the  gap,  we  first  note  that  the  calculated  offset  of  the 
valence-band  edge  turns  out  to  be  negligible,  due  to  the 
strong  anion  character  of  the  valence-band  maximum  in  the 
spl  tight-binding  scheme.  The  change  in  the  band  gap  is  al¬ 
most  entirely  accounted  for  by  the  lowering  of  the  conduc¬ 
tion-band  minimum  as  the  concentration  of  Hg  increases. 
Therefore,  the  change  of  the  band  gap  with  composition  has 
little  effect  on  the  position  of  the  deep  levels  measured  from 
the  valence-band  edge.  The  correct  band  offset  of  0.35  eVu 
may  be  obtained  by  including  the  effects  of  the  d  electrons 
However,  the  corrections  to  the  energies  of  the  deep  levels 
(relative  to  the  valence-band  edge)  due  to  d electrons  should 
lie  within  the  accuracy  of  tight-binding  calculations  such  as 
ours.  These  calculations  are  most  reliable  for  characteriza¬ 
tion  of  the  levels  and  energy  trends,  and  these  are  not  expect¬ 
ed  to  be  affected  by  the  inclusion  of  d  states. 

The  calculated  changes  in  the  energy  of  the  deep  levels  as 
the  concentration  of  Hg  in  Hg,  _xCd^Te  is  increased  result 
then  from  the  substitution  of  Hg,  instead  of  Cd,  on  most  of 
the  cation  sites  in  the  vicinity  of  the  interstitial.  The  shifts  of 
the  energy  levels  as  we  go  from  CdTe  to  Hg-rich  HgCdTe  are 
not  uniform.  For  tetrahedral  interstitials,  they  are  biggest  for 
defects  with  cation  nearest  neighbors,  and/or  a  sizable 
amount  of  electronic  charge  on  cations.  The  two  ^-symme¬ 
try  levels  shown  for  Cd  and  Hg  interstitials  satisfy  both  of 
these  conditions,  and  thus  they  exhibit  the  strongest  shift. 
This  leads  to  the  crossing  of  the  levels  in  the  intermediate 
concentration  range. 

Considering  the  self-interstitial  levels  which  could  be  re¬ 
sponsible  for  the  midgap  tunneling  in  Hg0  8  Cd0  2Te,  we  find 
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that  the  A  ,-symmetry  level  of  Hg  with  Te  nearest  neighbors 
(with  energy  £  =  £,,+  0.04  eV  in  CdTe)  appears  to  be  a 
good  candidate.  This  is  in  agreement  with  suggestions  which 
have  been  made,  based  on  experiment,  that  Hg  interstitials 
may  be  a  source  of  the  midgap  tunneling  levels.2  However, 
the  accuracy  of  the  present  calculation  is  not  sufficient  to 
determine  whether  this  A{  level  really  lies  in  the  gap  for 
HgorCd0  2Te.  The  locahzed  state  associated  with  this  Hg 
level  has  a  strong  cation  chatacter,  as  more  than  70%  of  the 
electronic  charge  is  located  on  cations.  Experimental  efforts 
to  further  characterize  the  midgap  tunneling  centers  would 
be  very  useful  in  identifying  these  centers. 

Indium  doping  was  found  to  give  rise  to  two  deep  levels  in 
the  gap  region  for  CdTe.  However,  neither  of  them  is  expect¬ 
ed  to  act  as  an  acceptor  level  which  could  lead  to  signficant 
self-compensation  in  CdTe.  The  hexagonal  In  level  (at 
£  =  Ev  4*  0.6  eV  in  CdTe)  is  occupied  for  a  neutral  intersti¬ 
tial,  and  thus  it  cannot  serve  as  an  acceptor  level  The  7V 
symmetry  tetrahedral  level  (at  E  =  £,.  -f  0.5  eV)  due  to  In 
with  cation  neatest  neighbors  is  partially  empty,  and  in  prin¬ 
ciple  it  could  serve  as  an  acceptor  level.  However,  it  seems 
more  likely  for  In  to  occupy  the  tetrahedral  position  between 
anions,  as  it  does  in  substitutional  sites.  Thus,  we  believe  that 
tetrahedral  In  with  cation  nearest  neighbors  is  energetically 
less  favorable  and  hence  much  less  numerous  than  tetrahe¬ 
dral  In  with  anion  nearest  neighbors.  The  latter  interstitial 
was  found  to  have  only  shallow  donor  levels,  as  does  substi¬ 
tutional  In.  Therefore,  we  conclude  that  In  in  tetrahedral  or 
hexagonal  interstitial  positions  does  not  introduce  acceptor 
levels  which  could  account  for  the  electrical  inactivity  found 
in  MBE-grown  In-doped  CdTe.  The  formation  energy  cal- 
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culations  which  are  needed  for  the  final  confirmation  of  our 
assumption  concerning  the  most  favorable  tetrahedral  site 
are  now  in  progress.16 
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The  previous  simple  model  used  to  calculate  energy  trends  for  the  deep  levels  due  to  localized 
electronic  states  associated  with  interstitial  and  substitutional  defects  has  been  extended  to  obtain 
estimates  of  the  total  energy  for  interstitials  at  different  sites  in  the  lattice.  These  calculations  have 
been  applied  to  predict  the  preferred  tetrahedral  site  for  Hg,  Cd,  In,  and  Te  interstitials  in 
Hg,  _  xCdtTe,  and  to  investigate  a  possible  diffusion  path  for  these  interstitials. 


I.  INTRODUCTION 

Both  the  electronic  behavior  and  the  compositional  profile 
resulting  from  growth  or  processing  of  semiconductor  alloys 
can  be  strongly  influenced  by  the  presence  of  impurity  and 
self-interstitials.  Deep  levels  corresponding  to  localized  in¬ 
terstitial  states  may  lead  to  electrical  activity  for  interstitials 
with  a  sufficiently  favorable  formation  energy  to  be  numer¬ 
ous.  Even  when  other  defects  are  more  numerous,  intersti¬ 
tials  often  make  a  major  contribution  to  the  diffusion  pro¬ 
cesses  which  determine  dopant  profiles  and  to  chemical 
interdiffusion  in  superlattices  and  junctions.  Calculating  the 
total  energy  of  the  interstitials  at  different  positions,  and  the 
electronic  deep  levels  due  to  these  interstitials,  can  help 
broaden  our  understanding  of  the  contributions  of  intersti¬ 
tials  to  electronic  behavior  and  kinetic  processes  in  semicon¬ 
ductors. 

Calculations  of  the  localized  states  due  to  point  defects  in 
semiconductors  may  be  divided  into  two  main  categories. 
Green’s  function  calculations,  as  have  been  developed  for 
interstitial  defects  in  Si,1  account  fully  for  the  periodicity  of 
the  host  crystal.  Methods  in  which  a  localized  basis  is  sepa¬ 
rated  from  the  rest  of  the  crystal,  or  repeated  periodically, 
are  only  valid  when  it  can  be  shown  that  the  boundary  condi¬ 
tions  chosen  for  the  localized  basis  do  not  strongly  alter  the 
description  of  the  localized  states  from  what  would  be  ex¬ 
pected  for  an  isolated  defect  in  an  infinite  crystal.  However, 
these  methods  are  more  easily  applicable  to  problems  such  as 
defects  with  lattice  distortions  and  defects  in  alloys. 

We  have  previously  developed  a  simple  tight-binding 
model  for  interstitial  and  substitutional  defects  in  tetrahe- 
drally  bonded  semiconductors  with  significant  ionic  and  me¬ 
tallic  character.2  In  this  approach,  the  localized  states  are 
described  '.n  terms  of  a  small  localized  basis.3  For  substitu¬ 
tional  and  high-symmetry  (tetrahedral)  interstitials,  the  en¬ 
ergy  levels  and  wave  functions  of  the  localized  states  are 
obtained  analytically.  Energy  levels  and  wave  functions  have 
been  obtained  numerically  for  hexagonal  interstitials.  Al¬ 
though  the  calculated  deep  levels  cannot  be  viewed  as  quan¬ 
titatively  reliable,  the  general  energy-level  structure  ob¬ 
tained  from  these  simple  calculations  compares  well,  with  a 
few  exceptions,  with  self-consistent  Green's  function  and  ex¬ 
perimental  results  for  interstitials  and  substitutional  defects 


in  Si,  GaP,  GaAs,  and  CdTe.2  Correct  symmetries  and  occu¬ 
pancies  are  obtained,  and  known  chemical  trends  repro¬ 
duced.  In  this  work,  we  investigate  the  effects  of  charge 
transfer  on  the  deep-level  structure,  and  extend  these  calcu¬ 
lations  to  obtain  estimates  of  the  total  energies  for  intersti¬ 
tials  at  different  sites  in  the  lattice. 

While  experimental  evidence  indicates  that  interstitials 
may  be  correlated  with  the  midgap  tunneling  levels  in  Hg- 
rich  HgCdTe,4  reliable  identification  of  specific  defect  levels 
in  HgCdTe  has  been  a  very  difficult  problem.  Previous  re¬ 
sults  of  this  model  for  the  deep  levels  due  to  In  and  self¬ 
interstitials  in  Hg,  _  ,CdxTe  alloys  have  been  reported,  and 
implications  for  self-compensation  and  doping  behavior  dis¬ 
cussed,  assuming  that  the  preferred  tetrahedral  site  for  Hg 
and  In  interstitials  is  between  anions.2  Total  energies  are 
calculated  here  for  In  and  Hg  self-interstitials  at  different 
sites,  and  the  preferred  tetrahedral  site  for  each  interstitial  is 
determined. 

Experimental  evidence  also  indicates  that  interstitial  de¬ 
fects  play  a  major  role  in  cation  interdiffusion  at  CdTe/ 
HgTe  junctions.5  Developing  a  more  accurate  picture  of  in¬ 
terstitials  in  these  materials  is  therefore  as  important  for  un¬ 
derstanding  diffusion  mechanisms  as  it  is  for  understanding 
the  microscopic  structure  and  properties  of  defects  associat¬ 
ed  with  the  midgap  levels  and  other  observed  electronic  lev¬ 
els.  Results  of  the  total  energy  calculations  for  In  and  self¬ 
interstitials  are  used  to  investigate  whether  the 
low-electron-density  channel  which  passes  alternately 
through  tetrahedral  sites  between  cations  and  tetrahedral 
sites  between  anions  is  likely  to  be  an  important  diffusion 
path  for  these  interstitials  in  HgCdTe.  Possibilities  for 
charge-transfer-assisted  diffusion  are  discussed. 

II.  LOCAL  MATRIX  APPROACH 
A.  Basic  model 

In  this  approach,  the  localized  states  associated  with  the 
defect  are  described  in  terms  of  a  small  basis  of  atomic  orbi¬ 
tals,  localized  around  the  defect.  These  basis  states  are  de¬ 
coupled  from  the  rest  of  the  crystal  by  the  method  of  “soft 
separation,”  i.e.,  without  breaking  any  of  the  covalent  bonds 
(the  strongest  coupling  in  the  crystal).  The  interactions  of 
the  interstitial  with  the  rest  of  the  localized  basis  states  arc 
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treated  exactly,  within  the  two-center  tight-binding  approxi¬ 
mation,  while  the  weaker  interactions  which  couple  the  local 
basis  to  the  rest  of  the  crystal  are  systematically  neglected. 
The  effects  of  the  polar,  metallic,  and  covalent  coupling 
within  the  localized  basis  are  thus  included,  while  preserving 
the  simplicity  of  the  small-basis  description.2 

The  minimal  basis  for  the  localized  interstitial  states  in 
this  model  includes  the  atomic  orbitals  on  all  the  nearest 
neighbors  of  the  interstitial,  which  are  improperly  coordi¬ 
nated  due  to  the  presence  of  the  defect,  and  all  the  sp*  hybrids 
on  more  distant  atoms  which  are  involved  in  bonds  with 
hybrids  on  the  nearest-neighbor  atoms.  Adding  the  atomic 
orbitals  on  the  interstitial  atom  itself  gives  a  basis  of  36 
states,  composed  of  atomic  orbitals  on  15  atoms,  for  the  te¬ 
trahedral  interstitial,  and  a  basis  of  40  states,  composed  of 
atomic  orbitals  on  19  atoms,  for  the  hexagonal  interstitial. 
The  analogous  basis  for  substitutional  defects  contains  32 
states,  composed  of  atomic  orbitals  on  17  atoms.  The  local 
structure  around  each  of  these  defects,  including  all  the 
atoms  which  contribute  states  to  the  localized  basis,  is  shown 
in  Fig.  1. 

The  minimal  basis  is  cut  off  from  the  rest  of  the  crystal  by 
neglecting  the  metallic  coupling  (which  is  proportional  to 
the  promotion  energy,  Ep  —  Es)  on  the  atoms  bonded  to  the 
nearest  neighbors,  and  all  the  two-center  interactions 
between  hybrids  on  the  nearest-neighbor  atoms  and  hybrids 
on  further  atoms,  except  the  ionic  and  covalent  coupling 
between  the  hybrids  on  the  nearest-neighbor  and  further 
atoms  which  are  covalently  bonded  to  each  other.  This 
method  of  separation  does  not  introduce  surface  states  due 
to  dangling  bonds  or  bonds  of  a  different  type,  such  as  the 
host-hydrogen  bonds  which  occur  when  the  bonds  at  the 
surface  of  the  local  cluster  are  terminated  with  hydrogen. 


B.  Tight-binding  parameters  and  approximations 

The  deep-level  energies  and  wave  functions  given  by  the 
formulas  derived  for  substitutional  defects  and  tetrahedral 
interstitials,2  and  found  numerically  for  hexagonal  intersti¬ 
tials,  depend  on  the  approximations  and  the  tight-binding 
parameters  used  to  define  the  matrix  elements  of  the  Hamil¬ 
tonian.  Modifications  in  the  set  of  parameters  used  will  cause 
some  change  in  the  calculated  levels.  However,  recent  work 
on  structural  total  energies  for  bulk  semiconductors6  has 
shown  that  the  results  of  such  tight-binding  total  energy  cal¬ 
culations  are  quite  insensitive  to  modifications,  provided 
that  the  parameters  reflect  the  chemical  trends  in  the  atomic 
data.  In  this  work,  we  follow  the  ideas  of  Harrison,7  and 
assume  that  the  on-site  matrix  elements  are  most  sensitive  to 
the  chemical  character  of  the  atom  occupying  the  site,  while 
off-site  matrix  elements  are  most  sensitive  to  the  interatomic 
distance. 

The  on-site  energies  are  determined  using  atomic  orbital 
energies  obtained  from  experiment,  supplemented  by  Har- 
tree-Fock  results  where  necessary.  The  first  ionization  ener¬ 
gies  taken  from  x-ray  data*  are  used  to  determine  the p-orbi- 
tal  energy  Ep  for  atoms  with  an  occupied  p  valence  shell,  and 
to  determine  the  j-orbital  energy  E \  for  atoms  with  only  s 
valence  electrons  (i.e.,  Hg,  Cd).  For  atoms  with  only  s  va- 


O  defect  atom 
O  cation 
#  anion 

Fig.  I.  Atomic  structure  of  the  (a)  substitutional  defect.  < b )  tetrahedral 
interstitial,  and  (c)  hexagonal  interstitial. 


lence  electrons,  the  average  excitation  energies  are  also  avail¬ 
able  from  the  x-ray  data  (and  are  not  available  from  Har- 
tree-Fock  calculations),  and  these  arc  used  to  determine  the 
promotion  energies.  Ep  -  Ey .  The  average  excitation  ener¬ 
gies  are  generally  not  available  from  the  experimental  data 
for  atoms  with  occupied p  valence  shells.  Therefore,  for  these 
atoms  the  promotion  energy  is  taken  from  Hartree-Fock 
results/'  The  resulting  energies  for  the  atoms  considered  in 
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this  work  are  Ex  =  -  10.43,  -  8.99,  -  18.28.  and  -  10.56 
cV  and  Ep  =  -  5.00,  -4.77,  -9.01.  and  -  5.78  eV  for 
Hg,  Cd,  Te,  and  In,  respectively. 

We  use  the  two-center  approximation  for  the  off-site  ma¬ 
trix  elements,  and  assume  that  these  matrix  elements  depend 
only  on  the  bond  length  d.  Following  Harrison,7  we  assume 
that  the  matrix  elements  describing  the  interaction  between 
nearest-neighbor  atoms  have  the  form  Vu  mfr/md % 

with  the  tight-binding  parameters  ?/n<,  =  -  1.32, 

—  1  42,  i]rpi1  —  2,22.  and  =  -  0.63,  where  m  is  the 
free-electron  mass.  Although  these  coefficients  were  devel¬ 
oped  by  Harrison  to  describe  the  interactions  between  atoms 
in  perfect  tetrahedrally  bonded  crystals,  we  assume  that  the 
interaction  between  the  interstitial  and  the  host  atoms  has 
the  same  dependence  on  d ,  and  may  be  described  by  the  same 
set  of  parameters. 

The  tight  binding  band  structure  obtained  for  perfect 
crystals  has  ofte  i  been  improved  by  adding  an  excited  s* 
orbital7,10  or  small  next-nearcst-neighbor  interactions/* 
which  could  be  included  in  the  present  small-basis  model. 
We  would  like  to  treat  the  chemical  trends  and  the  bonding 
for  host  and  interstitial  atoms  on  the  same  footing.  However, 
the  justification  for  including  an  orbital  or  second-near¬ 
est-neighbor  interactions  for  the  interstitial,  using  param¬ 
eters  fit  to  the  bulk  band  structure,  is  not  obvious.  Therefore, 
we  use  only  nearest-neighbor  interactions  in  the  present 
work. 


Ml.  CHARGE  TRANSFER  EFFECTS 

We  note  that  the  effective  intra-atomic  repulsion  for 
charge  transfer  between  cations  and  anions  in  the  perfect 
crystal  is  close  to  zero  for  the  materials  considered,  as  dis¬ 
cussed  in  detail  by  Harrison. n  We  therefore  neglect  the  ef¬ 
fect  of  charge  transfer  on  the  on-site  matrix  elements  for  the 
atoms  in  the  perfect  crystal.  When  the  interstitial  is  intro¬ 
duced,  additional  charge  transfer  between  the  interstitial 
and  the  host  atoms  occurs.  The  effects  of  this  charge  transfer 
to  or  from  the  interstitial  should  be  greatest  on  the  intersti¬ 
tial  itself. 

We  estimate  the  effects  of  this  charge  transfer  by  introduc¬ 
ing  an  effective  intra-atomic  Coulomb  repulsion  between 
electrons  on  the  intersitital.  The  Coulomb  shift  resulting 
from  removal  of  an  electron  from  the  interstitial  atom  to  an 
infinite  distance,  —  £/,  is  obtained  from  the  difference  of  the 
first  and  second  ionization  energies  for  this  atom,  taken  from 
x-ray  data.8  In  the  case  of  In,  where  the  first  electron  is  re¬ 
moved  from  the p  shell,  and  the  second  from  the  s  snell.  U  is 
obtained  from  x-ray  data,8  as  is  done  for  all  other  interstitial 
atoms.  However,  the  difference  between  the  first  and  second 
ionization  energies  for  In  is  not  equal  to  —  U ,  but  to 
Es—  Ep  —  U.  In  these  calculations,  charge  transferred  from 
the  neutral  interstitial  is  not  removed  to  an  infinite  distance, 
but  is  tranferred  to  the  nearest  neighbors  of  the  interstitial,  at 
a  distance  of  d(i,  or  to  atoms  bonded  to  the  nearest  neighbors, 
which  are  not  much  further  away.  Therefore,  the  effective 
Coulomb  shift  resulting  from  transfer  of  an  electron  between 
the  interstitial  and  nearby  host  atoms  is  given  by  Uctt 
=  U-  \4AeV(A/d0). 


After  the  Ham/  onian  matrix  is  diagonalized  to  obtain 
the  deep-level  enemies  and  the  wave  functions  of  the  local¬ 
ized  states,  the  shift  in  the  on-site  energy  for  the  interstitial 
due  to  charge  transfer  to  the  surrounding  atoms  is  calculat¬ 
ed,  and  the  new  Hamiltonian  matrix,  including  this  energy 
shift,  is  diagonalized.  This  procedure  is  repeated  until  a  self- 
consistent  solution  is  reached.  The  resi  lung  charge  on  the 
interstitial  itself,  before  the  defect  becor  es  ionized,  is  shown 
in  Table  I  for  various  interstitials  at  different  sites  in  Of  e 
and  Hg-rich  HgCdTe.  Interstititais  in  the  tetrahedral  posi¬ 
tion  between  cations  (T1 ),  where  the  electron  density  is  low, 
tend  to  lose  some  electronic  charge  to  the  surrounding 
atoms.  Interstitials  in  the  tetrahedral  position  between  an¬ 
ions  (T2).  where  the  electron  density  is  higher,  tend  to  gam 
some  electonic  charge  from  their  neighbors.  These  charge 
transfers  are  on  the  order  of  25%  of  an  electron.  Interstitials 
in  the  hexagonal  position,  which  have  an  equal  number  of 
cation  and  anion  nearest  neighbors,  tend  to  have  much 
smaller  charge  transfers.  The  preferred  ionization  states  for 
these  interstitials,  assuming  a  Fermi  level  near  the  middle  of 
the  gap.  are  also  listed  in  Table  I.  Interstitials  in  the  T1  posi¬ 
tion.  in  the  region  of  low  electron  density,  tend  to  remain 
neutral  or  become  negatively  charged,  while  interstitials  m 
the  T2  and  hexagonal  positions  tend  to  become  positively 
charged. 

The  shifts  in  the  deep-level  energies  for  neutral  intersti¬ 
tials  due  to  charge  transfer  are  on  the  order  of  0.4  eV  for 
tetrahedral  interstitials,  and  0.05  eV  for  hexagonal  intersti¬ 
tials.  The  effects  of  charge  transfer  do  not  change  the  qualita¬ 
tive  conclusions  reported  earlier.2  There  is  still  a  level  of  A , 
symmetry  in  the  vicinity  of  the  gap,  due  to  Hg  interstitials  in 
the  T2  position.  This  level  is  of  interest,  since  it  has  been 
suggested4  that  the  midgap  tunneling  levels  in  Hg-rich 
HgCdTe  are  associated  with  Hg  interstitials.  If  the  T2  posi¬ 
tion  for  Hg  interstitials  is  preferred,  it  is  certainly  possible 
that  the  midgap  tunneling  levels  of  Hg-rich  HgCdTe  arc  due 


Tabix  1.  Charge  on  interstitials. 


Material 

Interstitial  Position 

Preferred 

charge 

state 

Net  charge 
on  interstitial 
before  ionization 
(units of  c  ) 

CdTe 

Cd  Tl- 

Cd" 

-r  0.3 

llh 

Cd*- 

0. 

T2 

Cd*  ‘ 

-  0.3 

in  II 

In" 

-0.4 

H 

In’  ‘ 

-0.1 

T2 

In  ’ 

-  0.2 

Te  T1 

Te 

-  0.3 

H 

Te'  ’ 

0. 

T2 

Tc  ‘ 1 

-  0.1 

Hg  rich 

Hg  Tl 

Hg" 

-0.2 

HgCdTe 

H 

Hg  - 

-  0.3 

T2 

Hg 

-0.2 

’'1 1  s tetrahedral  position  between  cations. 
*‘Hs hexagonal  position. 

T2  =  tetrahedral  position  between  anions. 
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to  T2  Hg  interstitials,  pei  haps  associated  with  another  defect 
or  complex  which  pulls  this  level  into  the  gap. 

The  deep-level  energies  for  the  In  interstitial  in  the  T2 
position  are  not  strongly  affected  by  charge  transfer,  and  any 
shifts  which  do  result  from  charge  transfer  in  this  case  only 
tend  to  make  the  T2  In  interstitial  an  even  stronger  donor. 
Therefore,  as  previously  discussed,2  if  the  T2  position  is  pre¬ 
ferred  for  the  In  interstitial,  we  do  not  expect  self-compensa¬ 
tion  in  In-doped  HgCdTe  due  to  isolated  In  interstitials.  The 
electrical  inactivity  which  is  often  observed  in  In-doped 
CdTe  films  grown  by  MBE  without  laser  illumination12 
must  result  from  the  failure  of  the  In  to  occupy  isolated  sub¬ 
stitutional  sites,  or  compensation  of  the  shallow  substitu¬ 
tional  donors  by  some  other  defects. 


IV.  OVERLAP  AND  TOTAL  ENERGY 


In  order  to  calculate  total  energies,  the  short-range  repul¬ 
sion  arising  from  overlap  of  the  electron  wave  functions  cen- 
ered  on  neighboring  atoms  must  be  included.  The  overlap 
was  assumed  to  be  proportional  to  the  hopping  matrix  ele¬ 
ment  divided  by  the  average  on-site  energy,  and  to  have  the 
form  assumed  by  Majewski  and  Vogf: 


Sll'm  —  ’ 


1 

2 

£ah 

.  £/a  £/*u  . 

md *' 

where  is  the  atomic  orbital  energy  for  orbital  /  on  atom 
A,  and  r„  is  an  average  orbital  radius  determined  from  the 
Coulomb  repulsion  energy  of  two  electrons  on  the  same 
atom  (C/A  for  atom  A): 


*AB  is  a  scaling  constant  close  to  1,  which  depends  on  the 
row  of  the  Periodic  Table  containing  A  and  B.  For  A  and  B 
from  different  rows,  A’aij  is  given  by  the  geometric  average. 
Since  we  use  values  for  the  atomic-orbital  energies  and  for  V 
determined  from  experiment,  and  Harrison’s  tight-binding 
parameters,  while  Majewski  and  Vogl  use  Hartree- 
Fock  results  for  the  orbital  energies,  values  obtained  from 
Harrison’s  universal  parameter  theory  for  L\  and  a  slightly 
different  set  of  tight-binding  parameters,  we  have  deter¬ 
mined  the  values  of  A*  appropriate  for  our  set  of  parameters 
by  fitting  to  the  equilibrium  lattice  constants  for  Si,  C  and 
Sn.  This  procedure  gives  us  values  for  K  which  are  j  of  the 
values  of  Majewski  and  Vogl. 

If  we  approximately  orthogonalize  the  localized  basis 
functions,  taking  into  account  their  o\  erlap  to  first  order,  w  e 
will  get  corrections  to  both  the  off-site  and  the  on-site  matrix 
elements  of  the  Hamiltonian.  However,  since  the  off-site  ma¬ 
trix  elements  were  constructed  empirically,  they  should  al¬ 
ready  include  the  effects  of  overlap.  The  shifts  in  atomie- 
crbial  energies  resulting  from  overlap  are  typically  the  same, 
to  within  ±  30%,  for  all  the  orbitals  in  the  minimal  basis  for 
a  particular  interstitial  defect.  Therefore,  we  have  taken  the 
effect  of  overlap  to  be  a  constant  energy  shift  equal  to  the 
average  energy  shift  for  the  occupie  i  srbitals  in  the  minimal 
basis  for  each  defect. 

The  total  energy  for  a  semiconductor,  with  or  without  an 
interstitial  defect,  is  given  by 
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where  the  sum  of  electron  energies  €n  is  over  the  occupied 
one-electron  states,  the  electron-electron  interactions 
(which  have  been  included  twice  in  the  first  term)  are  sub¬ 
tracted  out,  and  the  core-core  interactions  are  added.  If  the 
material  is  composed  of  spherically  symmetric,  neutral 
atoms  with  no  overlap  of  charge  densities,  the  last  two  terms 
cancel. 

We  may  assume  that  these  two  terms  cancel  for  the  mate¬ 
rials  we  are  considering,  in  the  absence  of  charge  transfer. 
When  electrons  are  transferred  and  the  atoms  become 
charged,  the  electron-electron  interaction  energy  changes. 
In  setting  up  the  Hamiltonian  to  solve  for  the  one-electron 
energy  levels  €n ,  we  have  assumed  that  the  only  change  in  the 
electron-electron  interaction  energy  which  occurs  as  a  result 
of  charge  transfer  is  due  to  the  effective  Coulomb  repulsion 
which  opposes  the  transfer  of  more  electrons  to  the  intersti¬ 
tial  from  its  immediate  vicinity.  When  the  electronic  charge 
on  the  interstitial  changes  from  the  number  of  valence  elec¬ 
trons  on  the  neutral  atom  Z,.  to  (?,.  the  resulting  change  in 
the  electron-electron  interaction  energy  is 
UKtr (<2;  —  Z2)/2,  which  must  be  subtracted  from  the  sum 
of  the  energies  of  the  occupied  one-electron  states  to  give  the 
total  energy. 

Effects  of  relaxation  and  lattice  distortion  have  been  ne¬ 
glected.  If  relaxation  of  the  lattice  is  allowed,  the  total  ener¬ 
gies  may  be  modified.  However,  lattice  relaxation  is  expect¬ 
ed  to  be  less  important  for  tetrahedral  interstitials  than  for 
interstitials  at  sites  with  lower  symmetry  and  less  favorable 
nearest-neighbor  distances.  Therefore,  inclusion  of  these  ef¬ 
fects  shold  not  change  the  trends  discussed  in  the  present 
work. 

We  report  here  the  total  energies  of  the  localized  basis, 
containing  a  neutral  interstitial  defect,  relative  to  the  total 
energy  of  the  localized  basis  with  no  interstitial,  plus  the 
isolated  interstitial  atom  at  an  infinite  distance.  The  iota! 
energies  obtained  from  this  model  are  not  expected  to  be 
quantitatively  accurate,  but  provide  estimates  useful  for  in¬ 
dicating  the  preferred  position  of  interstitials,  and  studying 
trends  in  the  barriers  to  interstitial  diffusion  along  different 
paths.  Using  the  same  approximations  for  the  total  energy  as 
described  above,  we  obtain  calculated  values  for  the  cohesive 
energy  of  1.4eV/bond  for  both  CdTe  and  HgTc,  which  may 
be  compared  with  the  experimental  values  of  1.1  eV/bond 
(CdTe)  and  0.8  eV/bond  (HgTe).  The  bonding  energy  is 
overestimated  in  HgTe.  due  to  neglect  of  the  effect  of  the  d 
states  on  the  bonding. 


V.  PREFERRED  INTERSTITIAL  SITES 

Total  energy  differences  for  selected  interstitials  at  differ¬ 
ent  sites  in  the  latticcarc  given  in  Table  II.  ,  —  £',\  is  the 
difference  of  the  total  energies  for  the  neutral  interstitial  in 
the  T1  and  T2  tetrahedral  positions,  requiring  all  bound 
electrons  or  holes  to  be  in  localized  deep  levels  oi  .  csonances. 
E*\  i  —  £',*:  is  the  difference  of  the  total  energies  for  the  neu¬ 
tral  interstitial  in  theTl  and  T2  tetrahedral  positions,  allow¬ 
ing  for  the  formation  of  shallow  hydrogenie  u  -r  or  accep* 
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Table  II.  Interstitial  total  energies/ 


po  pO 

£  ti  *" 

|T°  _  pO 

£T1  —  Ejy 

Preferred 

Material 

Interstitial 

(eV) 

(eV) 

(eV) 

transition 

CdTe 

Cd 

+  6.4 

+  7.7 

+  10.9  -2 fi 

Cti'ri  ~CdTy 

In 

+  1.3 

+  3.2 

+  4.8  -  ft 

In’/i  Inf: 

Te 

-  1.8 

(  —  3.5)h 

(  +2.9-6 fi)* 

Tcrf-Tef/ 

Hg  rich 

HgCdTe 

Hg 

-  1.8 

+  1.2 

+  \A-2fi 

Hg?i-HgTy 

4  T1  =  tetrahedral  positton  between  cations.  T2s  tetrahedral  position  between  anions.  total  energy  for  neutral  interstitial  with  all  electrons  in  localized 
deep  levels  or  resonances.  £u  -  total  energy  for  neutral  interstitial  allowing  shallow  donor/acceptor  levels.  £=  total  energy  for  preferred  charge  state. 

=  Fermi  level,  measured  from  the  valence-band  edge. 
h  Not  accurate  due  to  large  charge  transfer. 


tor  levels,  so  that  bound  electrons  or  holes  may  occupy  these 
hydrogenic  states,  with  energies  near  the  band  edge,  rather 
than  deep  localized  states  with  energies  outside  the  gap. 
Eri  —  Eri  is  the  difference  of  the  total  energies  for  intersti¬ 
tials  in  the  preferred  ionization  state,  listed  in  Table  I,  forTl 
and  T2  tetrahedral  positions.  The  deep-level  structure  ob¬ 
tained  in  these  self-consistent  calculations  for  neutral  inter¬ 
stitials  does  not  provide  a  good  representation  for  highly 
ionized  interstitials,  due  to  the  large  charge  transfer.  For  this 
reason,  we  have  listed  the  total  energy  differences  between 
Te  interstitials  in  TI  and  T2  positions  in  parentheses,  when 
the  many  holes  and  electrons  in  deep  resonant  states  have 
been  allowed  to  fall  to  shallow  hydrogenic  bound  states,  or 
the  interstitials  have  ionized  to  their  preferred  charge  state. 

As  shown  in  Table  I,  the  calculated  total  energy  for  neu¬ 
tral  indium  interstitials  is  higher  in  the  tetrahedral  position 
between  cations  (Tl)  than  in  the  tetrahedral  position 
between  anions  (T2)  by  -3  eV  in  CdTe.  The  preferred 
charge  state  of  these  interstitials  depends  on  the  Fermi  level 
/i,  but  for  a  Fermi  level  of  -  1  eV,  measured  from  the  va¬ 
lence-band  edge,  the  preferred  charge  state  is  neutral,  for  the 
Tl  In  interstitial,  and  +  1  for  theT2  In  interstitial.  Since  the 
T2  interstitial  has  an  even  lower  energy  when  it  ionizes  to 
In  * ,  the  total  energy  difference  between  Tl  and  T2  In  inter¬ 
stitials  in  their  preferred  charge  states  is  even  greater  than  for 
neutral  interstitials.  These  results  support  our  earlier  as¬ 
sumptions2  that  the  preferred  tetrahedral  site  for  In  intersti¬ 
tials  is  T2.  Therefore,  although  Tl  In  interstitials  have  a 
partially  filled  deep  level  in  the  gap  and  could  potentially 
serve  as  acceptors,2  these  interstitials  will  not  be  numerous 
enough  to  account  for  the  electrical  compensation  which  is 
seen  in  non-laser-assisted  MBE-grown  films  of  CdTe.12 

Vydyanath  has  looked  at  the  effect  of  Hg  overpressures  on 
the  midgap  tunneling  centers  in  Hg-rich  HgCdTe,  and  sug¬ 
gests  these  centers  may  be  due  to  Hg  interstitials.4  We  have 
previously  found  that  the  T2  Hg  interstitial  has  an  A ,  level  in 
the  gap  region  2  In  this  work,  we  find  that  the  total  energy  is 
lower  for  Hg  interstitials  in  T2  than  in  Tl  sites.  The  energy 
difference  is  -  I  eV  for  neutral  interstitials,  and  slightly 
more  when  the  T2  interstitial  ionizes  to  Hg*2.  The  total 
energy  differences  listed  for  Hg-rich  HgCdTe  in  Table  II  are 
for  material  with  a  band  gap  of  0. 1  eV,  and  a  Fermi  level  ft 
-0.05  eV.  According  to  these  results,  tetrahedral  Hg  inter- 
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stitials  will  tend  to  sit  in  T2  sites,  and  contribute  Ax  elec¬ 
tronic  levels  in  the  vicinity  of  the  gap. 

The  total  energy  differences  listed  in  Table  II  for  Cd  and 
Te  interstitials  in  CdTe  follow  the  same  general  trend  as  was 
verified  above  for  Hg  and  In  interstitials,  elements  which  sit 
primarily  on  cation  substitutional  sites  will  prefer  the  tetra¬ 
hedral  interstitial  position  on  the  alternate  cation  sublattice, 
and  elements  which  sit  primarily  on  anion  lattice  sites  will 
prefer  the  tetrahedral  position  on  the  alternate  anion  sublat¬ 
tice.  We  find  that  the  T2  position  (between  anions)  is  pre¬ 
ferred  for  tetrahedral  Cd  interstitials,  while  the  Tl  position 
(between  cations)  is  preferred  for  tetrahedral  Te  intersti¬ 
tials. 

VI.  DIFFUSION 

Total  energy  differences  for  interstitials  at  different  sites 
can  provide  useful  estimates  of  diffusion  barriers,  and  help 
identify  probable  paths  for  interstitial  diffusion.  If  one  of  the 
tetrahedral  sites  is  the  preferred  interstitial  position,  the  path 
through  the  channels  of  low  electron  density,  passing  alter¬ 
nately  through  Tl  and  T2  tetrahedral  sites,  may  be  an  im¬ 
portant  avenue  for  diffusion.  We  may  begin  to  investigate 
diffusion  barriers  along  this  path  by  examining  the  energy 
differences  for  interstitials  in  Tl  and  T2  sites. 

Experimental  studies  in  both  CdTe  and  Hgo^Cd^Te 
suggest  that  interstitials  play  an  important  role  in  cation  dif¬ 
fusion,  and  that  Cd  interstitials  in  particular  play  an  impor¬ 
tant  role  in  interdiffusion  at  HgTe/CdTe  junctions/  The 
calculated  total  energy  difference  for  neutral  Cd  interstitials 
in  Tl  and  T2  positions  in  CdTe,  shown  in  Table  II,  is  —  8  eV, 
much  higher  than  experimental  activation  energies/  Simple 
interstitial  diffusion  through  the  low-electron-density  chan¬ 
nel  which  passes  alternately  through  T2  and  Tl  sites  there¬ 
fore  seems  unlikely  for  Cd. 

Calculated  total  energy  differences  for  neutral  interstitials 
in  Tl  and  T2  positions  are  —  3  eV  for  In  in  CdTe  and  —  1  eV 
for  Hg  in  Hg-rich  HgCdTe.  Although  the  energy  difference 
is  not  so  great  for  In  interstitials  as  for  Cd  interstitials  in 
CdTe,  the  open  channel  between  Tl  and  T2  tetrahedral  posi¬ 
tions  also  appears  to  be  unlikely  for  In  diffusion.  However, 
this  channel  appears  to  be  a  good  candidate  for  diffusion  oi 
Hg  in  Hg-rich  HgCdTe,  provided  that  the  energy  required  tu 
pass  through  the  hexagonal  site  between  tetrahedral  sites  is 
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not  excessive*  A  reliable  estimate  of  the  energy  required  for 
Hg  to  pass  through  the  hexagonal  site  may  require  the  inclu¬ 
sion  of  relaxation  outward  of  the  nearest  neighbors  of  the 
interstitial,  which  is  beyond  the  scope  of  these  small-basis 
calculations. 

Diffusion  behavior  may  be  quite  different  for  interstitials 
which  are  not  neutral.  As  shown  in  Table  I.  the  preferred 
charge  state  for  all  the  interstitials  listed  changes  radically  at 
different  sites.  Cation  interstitials  in  HgCdTe  tend  to  be  neu¬ 
tral  in  T1  sites,  and  4-  2  in  T2  sites  and  hexagonal  sites.  In 
interstitials  tend  to  be  neutral  in  T1  sites  and  4-  1  to  4-  3  in 
hexagonal  and  T2  sites.  Anion  interstitials  tend  to  be  -  2  in 
T1  sites  and  4-  2  to  -f  4  in  T2  and  hexagonal  sites.  This 
suggests  that  the  alternating  capture  and  release  of  charge 
carriers  by  interstitials  may  be  important  in  lowering  energy 
barriers  for  interstitial  diffusion  in  these  materials.  If  inter¬ 
stitials  are  allowed  to  change  to  their  preferred  charge  state 
at  each  site,  the  diffusion  barriers  become  dependent  on  the 
Fermi  level.  However,  the  conclusions  stated  above  are  not 
changed:  the  open  channel  between  interstitial  sites  appears 
to  be  a  favorable  diffusion  path  for  Hg,  less  favorable  for  In, 
and  quite  unfavorable  for  Cd. 

VII.  CONCLUSION 

The  small-basis  approach  previously  de\  eloped  to  study 
trends  in  the  deep-level  energies  for  interstitial  and  substitu¬ 
tional  defects  has  been  extended  to  include  the  effects  of 
charge  transfer  and  overlap.  Estimates  of  the  total  energies 
of  interstitial  defects  in  the  tetrahedral  position  between  ca¬ 
tions  and  the  tetrahedral  psotion  between  amons  are  calcu¬ 
lated  within  this  model.  The  preferred  tetrahedral  site  is 
found  to  be  between  cations  for  Te,  and  between  anions  for 
Cd,  Hg,  and  In.  Therefore,  as  previously  suggested. :  the  po¬ 
tential  acceptor  level  calculated  within  this  model  for  tetra¬ 


hedral  In  interstitials  between  cations  should  not  be  impor¬ 
tant  in  producing  the  electrical  compensation  often  observed 
in  In-doped  CdTe.12 

Diffusion  along  the  open  channel  between  tetrahedral 
sites  is  found  to  be  unfavorable  for  Cd  and  probably  unfavor¬ 
able  for  In  in  CdTe,  but  may  be  important  for  Hg  in  Hg-nch 
HgCdTe.  Charge  carrier  capture  is  predicted  to  lower  diffu¬ 
sion  barriers  for  all  these  interstitials. 
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APPENDIX  C  C-i 

Point  defects  with  lattice  distortion  in  CdTe  and  HgCdTe 

J.T.  Schick  and  C.  G.  Morgan-Pond 
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( Received  4  October  1989;  accepted  29  October  1 989 ) 

Results  of  self-consistent,  tight-binding  supercell  calculations  of  electronic  defect  levels  and 
relaxation  about  defects  for  As.  Sb,  and  Te  interstitials  in  CdTe  and  for  the  Hg  vacancy  in  HgTe 
are  presented.  Trends  in  the  localized  defect  levels  and  the  relaxation  of  the  lattice  seen  for  these 
defects  and  for  Cd,  Hg,  and  In  interstitials,  from  previous  calculations.'  are  discussed.  Te 
interstitials  are  found  to  have  the  largest  lattice  relaxation  of  these  interstitials  in  their  preferred 
tetrahedral  positions.  This  may  help  explain  the  tendency  for  excess  Te  to  condense  into  second- 
phase  inclusions.  It  is  also  suggested  that  a  strain-mediated  interaction  between  Te  interstitials 
and  cation  vacancies  may  promote  the  creation  of  Te  antisites  in  Te-rich  HgCdTe. 


I.  INTRODUCTION 

A  great  deal  of  attention  has  been  focused  on  the  formation 
of  defects  and  their  effects  on  semiconductor  devices.  De¬ 
fects  and  defect  complexes  can  fundamentally  affect  the  elec¬ 
tronic  properties  of  the  crystal.  In  addition,  diffusion  of  the 
defects  in  the  materia)  during  growth  has  a  profound  effect 
on  the  uniformity  of  the  crystal.  Even  at  room  temperature, 
defect  diffusion  may  play  a  significant  role  in  the  stability  of 
a  crystal  over  time. 

Technologically  important  defects  may  be  in  the  form  of 
simple  point  defects,  such  as  interstitials  or  vacancies,  or 
may  be  more  extensiye.  Lattice  relaxation  occurs  for  simple 
point  defects  as  yvell  as  for  more  extended  defects,  yvhich 
may  involve  substantial  lattice  distortion.  This  relaxation 
may  change  the  symmetry  of  the  defect  site,  and  generally 
does  change  the  charge  distribution  in  the  vicinity  of  the 
defect,  modifying  both  the  defect  level  energies  and  the  cor¬ 
responding  localized  defect  states.  The  resulting  strains  in 
the  lattice  may  affect  the  attraction  and  binding  of  point 
defects  into  a  complex,  or  to  extended  defects. 

Both  experimental  evidence2  and  theoretical  estimates  of 
diffusion  barriers'  indicate  that  Hg  interstitials  may  diffuse 
easily  along  the  open  channels  betyveen  tetrahedral  sites  in 
Hg,  ,  Cd  4  Te,  while  Cd  interstitial  diffusion  requires  a  more 
complicated  process,  involving  interchanges  with  atoms  on 
lattice  sites.  Ultimately,  calculation  of  total  energies  for 
highly  distorted  local  configurations  may  be  necessary  in 
order  to  identify  fa  zombie  paths  for  Cd  diffusion. 

First  principles  techniques  have  often  been  used  to  calcu¬ 
late  the  properties  of  simple  defects  yvith  only  a  feyv  degrees 
of  freedom.  For  example,  pseudopotential  density-func¬ 
tional  approaches4  and  self-consistent  Green's  function 
methods’  have  been  applied  yvith  reasonable  success  for  te- 
tranedrai  and  hexagonal  self-interstitials  m  silicon.  While 
these  methods  produce  the  most  reliable  results,  they  require 
copious  amounts  of  computer  time,  yvhich  is  undesirable 
when  the  number  of  problems  to  study  is  large,  or  yvhen 
many  degrees  of  freedom  must  be  included  in  order  to  find 
the  preferred  configuration  of  the  lattice  surrounding  a  de¬ 
fect.  Tight-binding  theory  offers  an  attractive  alternative  to 
first  principles  methods  because  of  the  reduced  computa¬ 
tional  effort  Transferability  of  the  parameters  to  systems 
other  than  perfect  crystals  has  been  demonstrated  by  Chadi0 


in  his  work  on  surfaces  of  group  IV  and  III-V  semiconduc¬ 
tors,  and  by  Tomanek7  in  his  yvork  on  Si  clusters.  Compari¬ 
son  of  tight-binding  results  to  the  results  of  first  principles 
methods4,5  for  interstitials  has  been  encouraging. s 

We  have  recently  developed  a  supercell  approach  for  the 
calculation  of  total  energies  in  systems  yvith  point  defects, 
based  on  the  tight-binding  model  of  Majeyvski  and  Vogi." 
and  have  applied  this  technique  to  self-interstitials  in  Siv  and 
Cd.  Hg,  and  In  interstitials  in  CdTe  and  HgTe.1  Lattice  re¬ 
laxation  is  included,  and  charge  transfer  is  treated  in  a  fully 
self-consistent  manner.  We  report  here  the  results  of  our 
studies  of  the  properties  of  As,  Sb.  and  Te  interstitials  in 
CdTe  and  of  the  Hg  vacancy  in  HgTe.  and  expand  our  inves¬ 
tigation  of  trends  in  the  electronic  levels  and  relaxation  asso¬ 
ciated  yvith  cation  and  anion  self-interstitials  and  dopant  in¬ 
terstitials  in  Hg,  xCdxTe. 

II.  THEORY  AND  PARAMETERS 

The  supercells  used  in  this  yvork  have  a  cubic  symmetry, 
yvith  a  side  length  equal  to  tyvice  the  lattice  constant  of  the 
crystal.  There  are  64  atoms  in  the  bulk  crystal  supercell,  one 
more  for  interstitials,  and  one  less  for  vacancies.  Only  ^  and p 
electronic  states  are  considered  to  participate  in  the  bonding 
for  CdTe  and  HgTe.  In  this  lattice,  an  interstitial  may  be 
placed  at  a  tetrahedral  position,  yvith  four  anion  (or  cation ) 
nearest  neighbors  and  six  cation  (or  anion)  next  nearest 
neighbors  only  slightly  further  ayvay.  The  hexagonal  posi¬ 
tion,  which  is  at  the  narrowest  passage  betyveen  tyvo  adjacent 
tetrahedral  points,  is  surrounded  by  a  puckered  ring  consist¬ 
ing  of  three  anions  and  three  cations.  A  modified  form  of 
Majeyvski  and  VoglV  energy  formalism  is  used  to  predict 
the  most  favorable  configuration  of  the  neighboring  atoms. 
We  permit  a  breathing  mode  of  relaxation  of  the  nearest 
neighbor  atoms  radially  toyvard  or  ayvay  from  the  defect,  and 
a  change  in  puckering  of  the  nearest  neighbor  ring  lor  the 
hexagonal  interstitials.  At  this  time  yve  have  not  yet  included 
Jahn-Teller  distortions,  yvhich  are  expected  for  the  neutral 
states  of  the  cation  vacancy  and  Te  and  dopant  interstitials, 
where  partially  occupied  degenerate  states  occur.  Jahn-Tell¬ 
er  distortion  will  not  I'^cur  for  charge  states  with  completely 
filled  or  empty  defect  levels,  such  as  the  doubly  charged  ca¬ 
tion  vacancy  in  HgCdTe. 
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The  interactions  are  given  by  Harrison’s  inverse  square 
tight-binding  form10 

hrm  —  Virm  — »  U ) 

ma  “ 

where  d  is  the  separation  between  the  interacting  atoms,  m  is 
the  mass  of  the  electron  and  7jlrm  is  a  “universal”  parameter 
for  the  interaction  between  the  / and  / '  orbitals  in  a  bond  of  m 
symmetry.  We  impose  a  cutoff  of  1.39  times  the  nearest 
neighbor  distance  in  the  bulk  crystal,  which  corresponds  to 
including  only  first  neighbor  interactions  in  the  bulk  crystal 
and  with  the  hexagonal  interstitial,  and  first  and  second 
neighbor  interactions  with  the  tetrahedral  interstitial. 

The  on-site  energies  as  expressed  by  Majewski  and  Vogl9 
are  given  by 
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where  wAl  is  the  energy  of  the  A  orbital  of  the  isolated,  neu¬ 
tral  atom  on  site  i.  As  the  charge  is  redistributed  in  ionic 
crystals  and  crystals  with  defects,  corrections  to  the  one  elec¬ 
tron  potential  resulting  from  this  charge  transfer  are  includ¬ 
ed.  The  quantity  (Z,  —  £>,)  represents  the  electron  defi¬ 
ciency  on  the  atom  at  site  /*,  where  Z  is  the  valence  of  the 
neutral  atom  and  Q  is  the  number  of  electrons  at  the  site.  The 
second  term  is  the  intra-atomic  correction,  and  the  third 
term  is  the  interatomic  correction.  Ut  is  the  average  Cou¬ 
lomb  repulsion  between  two  electrons  on  the  same  atom.11 
The  modified  Coulomb  potential,  used  in  the  third  term  to 
account  for  the  fact  that  at  small  separations  the  distinction 
between  intra-atomic  and  interatomic  is  meaningless,  is  giv¬ 
en  by10 
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where  e  is  the  charge  of  the  electron  and  V  is  the  average 
Coulomb  repulsion  for  the  atoms  in  the  crystal.  This  poten¬ 
tial  goes  to  U  as  the  distance  between  neighboring  atomic 
sites  becomes  small,  while  for  large  separations  it  becomes 
Coulombic.  The  summations  are  carried  out  for  the  third 
term  using  Ewald’s  method.12  Since  the  on-site  energies  de¬ 
pend  explicitly  on  the  distribution  of  charge  in  the  crystal,  it 
is  necessary  to  iterate  until  a  self-consistent  solution  for  the 
eigenstates  of  the  system  is  found. 

The  overlap  correction  [h  provides  the  repulsion  neces¬ 
sary  to  prevent  the  crystal  from  collapsing.  Including  effects 
of  nonorthogonality  to  first  order  in  the  overlap  S  this  cor¬ 
rection  is  ft,  =  —  1(577, i  -r  //c, S) where  HxX  is  the  one 
electron  Hamiltonian  made  up  of  the  interactions  and  on¬ 
site  energies  described  above  °  The  overlap  between  orbital  / 
on  atom  A  and  orbital  / '  on  atom  B  is n 
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where  K  is  a  numerical  constant  close  to  one,  and  /*„  is  the 
average  radius  of  a  valence  state  determined  from  the  Cou¬ 
lomb  repulsions  U  on  atoms  A  and  B. 


The  total  energy  of  the  crystal  (per  supercell)  is  given  by 
the  sum  of  the  one  electron  energies  of  the  occupied  states, 
corrected  for  the  double  counting  of  the  electron-electron 
interactions,  plus  a  contribution  due  to  interactions  between 
the  rigid  cores 

occ 
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The  electron-electron  energy  may  be  separated  into  two 
parts 

i  cell 

=  (6) 
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The  second  summation  in  Eq.  (7)  is  the  same  as  that  in  Eq. 
(2),  and  therefore  is  evaluated  in  the  same  way.  The  remain¬ 
ing  sums  in  Eqs.  (6)  and  (7)  involve  only  the  atomic  sites 
within  the  supercell. 

Harrison’s  parameters10  77 „„  =  -  1.40,  =  3.24, 7/rp. 

=  —  0.81,  and  =  1.84  were  used  for  the  interactions. 
Separately  defined  interactions  with  more  distant  neighbors 
are  not  included  in  our  current  calculations.  These  interac¬ 
tions  are  usually  small,9  and  cannot  be  transferred  unambig¬ 
uously  to  calculations  involving  interstitials. 

The  bare  atomic  5-orbital  energies  and  occupied  /7-orbital 
energies  are  obtained  from  Hartree-Fock  calculations,’4 
with  the  exception  of  Hg,  for  which  we  use  Dirac-Fock  ener¬ 
gies  15  The  excited  /7-level  energies  are  obtained  from  experi¬ 
mental  spectra.1*  The  energies  are,  in  eV,  =  —7.20, 
-  8.93.  -  10.14.  -15.83,  -  18.66,  and  -19.06  and 
=  -  2.99.  -  3.50,  -  5.37,  -  9.10,  -  10.05,  and  -9.79 
for  Cd,  Hg,  In,  Sb.  As,  and  Te,  respectively.  The  values  of  K 
and  thee/"'  dependence  of  the  overlap  S  were  determined  by 
Majewski  and  Vogl  to  give  the  best  agreement  with  the  bond 
lengths,  bulk  moduli,  and  preferred  crystal  structures  for  a 
wide  variety  of  insulators  and  semiconductors/’  For  the  row 
including  Cd,  In.  Sb,  and  Te,  K  =  1.05,  for  the  As  row, 
K  =  1.16,  and  for  the  Hg  row,  K  —  0.93.  K  is  given  by  the 
geometric  average  for  overlaps  between  atoms  in  different 
rows.  Finally,  the  values  of  the  intra-atomic  Coulomb  repul¬ 
sions.  determined  by  Harrison  from  experiment.1 1  are  (in 
eV)  7.33,  6.95,  6.00*7.39,  8.31,  and  8.00  for  Hg,  Cd,  In,  Sb, 
As,  and  Te,  respectively. 

Ill*  INTERSTITIAL  AND  VACANCY  LOCALIZED 
STATES 

Although  numerical  predictions  of  deep  defect  level  ener¬ 
gies  obtained  using  tight-binding  methods  are  not  quantita¬ 
tively  reliable,  the  general  energy  level  structure  obtained 
from  simple  tight-binding  calculations  using  a  small-basis 
“local  matrix*'  method  compares  well,  in  most  cases,  with 
experimental  and  self-consistent  Green's  function  results  for 
interstitial  and  substitutional  defects  in  CdTe.  GaAs,  GaP. 
and  Si.  1  In  this  work,  in  addition  to  a  more  complete  treat¬ 
ment  of  overlap  and  charge  transfer,  we  include  the  coupling 
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of  the  local  matrix  describing  the  defect  to  the  rest  of  the 
infinite  matrix,  by  considering  a  periodic  array  of  defects.  To 
determine  the  effect  of  the  supercell  on  the  localized  states, 
the  dispersion  of  these  states  along  the  [  100]  and  [  1 1 1  ]  di¬ 
rections  was  calculated.  The  dispersion  is  on  the  order  of 
0.05  eV,  except  for  the  localized  defect  states  close  to  the 
band  edges,  which  often  interact  strongly  »vith  bulk  states. 
The  resulting  states  are  of  mixed  character,  less  localized  on 
the  interstitial  and  showing  some  of  the  larger  dispersion 
expected  for  bulk  states  at  the  band  edge.  Dispersions  for 
these  states  can  be  up  to  0.2  eV.  In  a  few  of  these  cases, 
involving  strong  interactions  between  bulk  states  and  local¬ 
ized  defect  states  of  similar  energy,  it  appears  that  the  local¬ 
ized  character  shifts  between  different  states  as  we  move 
away  from  the  T  point  toward  the  zone  edge.  It  therefore 
appears  that  the  uncertainty  in  deep  defect  level  energies 
resulting  from  the  use  of  the  supercell  approach  is  less  than 
the  uncertainty  due  to  problems  inherent  in  the  tight-bind¬ 
ing  model. 

Results  for  the  deep  level  structure  associated  with  the 
interstitials  of  different  cation-hke  atoms,  such  as  Cd,  Hg, 
and  In,  which  prefer  to  sit  on  cation  sites  as  substitutional 
impurities,  follow  a  basic  pattern.  Cd  interstitials  in  the  pre¬ 
ferred  tetrahedral  position  (with  Te  nearest  neighbors)  in 
CdTe  have  a  defect  level  of  A ,  symmetry  close  to  the  conduc¬ 
tion  band  edge.  There  is  a  similar  level  of  A  symmetry, 
roughly  resonant  with  the  conduction  band  edge,  for  the 
hexagonal  Cd  interstitial,  and  a  comparable  A ,  level  for  the 
tetrahedral  Cd  interstitial  with  cation  nearest  neighbors. 
The/1 ,  level  for  the  tetrahedral  Hg  interstitial  with  Te  near¬ 
est  neighbors  is  somewhat  lower  in  energy,  and  this  level 
moves  downward  as  the  composition  of  the  crystal  is  varied 
from  CdTe  to  HgTe.  For  both  Hg  and  Cd  interstitials,  thc/i , 
level  holds  the  last  two  electrons  for  the  neutral  interstitial. 
For  Hg  and  Cd  interstitials  in  the  preferred  tetrahedral  posi¬ 
tion,  relaxation  of  the  lattice  around  the  interstitial  allows 
this  level  to  become  more  localized  on  the  interstitial  itself, 
and  to  move  slightly  lower  in  the  gap.  The  At  level  of  the 
tetrahedral  In  interstitial  with  Te  nearest  neighbors  falls  still 
lower  in  the  gap,  and  this  interstitial  also  has  a  7\  resonance 
well  above  the  conduction  band  edge,  w hich  donates  an  elec¬ 
tron  into  the  conduction  band. 

For  interstitials  of  anion-like  atoms,  such  as  Te,  As,  and 
Sb,  which  prefer  to  sit  on  anion  sites  as  substitutional  impuri¬ 
ties  in  the  absence  of  strong  deviations  from  stoichiometry ,  a 
similar  trend  is  observed  In  the  positions  of  the  localized 
states,  paralleling  the  trend  in  the  atomic  levels  of  the  isolat¬ 
ed  interstitial  atoms.  For  the  Sb  interstitial  between  cations 
in  CdTe,  there  is  a  partially  occupied  localized  state  of  Tm 
symmetry  high  in  the  gap,  and  an  Ax  state  close  to  the  va¬ 
lence  band  edge.  The  7\  level  for  the  As  interstitial  between 
cations  is  slightly  lower,  and  the  localized  A ,  state  is  now  a 
resonance  below  the  valence  band  edge.  For  the  Te  tetrahe¬ 
dral  interstitial  between  cations,  the  partially  occupied  7\ 
state  in  the  upper  part  of  the  gap  and  the  A ,  resonance  just 
inside  the  valence  band  are  again  slightly  lower. 

The  vacancy  is  modeled  as  an  atom  with  an  s  level  at  3000 
e V  and  a p  level  at  1 2  000  eV,  which  is  substituted  for  one  of 
the  atoms  in  the  bulk  lattice.  Because  of  the  range  of  the 


atomic  energies  involved,  results  for  vacancy  total  energy 
calculations  suffer  from  an  increased  level  of  computational 
error  resulting  in  an  increase  in  the  number  of  iterations 
needed  to  obtain  self-consistency.  For  the  Hg  vacancy  in 
HgTe,  a  localized  state  of  T2  symmetry,  with  over  half  of  its 
density  on  the  first  neighbors  of  the  vacancy,  is  found  near 
the  valence  band  edge.  Slightly  lower  in  energy  is  a  state  of 
Ax  symmetry,  which  is  even  more  strongly  localized  on  the 
nearest  neighbors  of  the  vacancy.  Similar  localized  states  of 
Ax  and  7\  symmetry  are  expected  to  be  near  the  valence 
band  edge  in  Hg-rich  HgCdTe. 

IV.  RELAXATION 

When  the  lattice  surrounding  the  defect  is  allowed  to  re¬ 
lax,  the  nearest  neighbors  of  the  interstitial  move  outward, 
and  the  nearest  righbors  of  the  vacancy  move  inward,  re¬ 
ducing  the  size  ot  .he  hole  left  at  the  vacant  site.  The  on-site 
energies  on  the  interstitial  atom,  and  often  on  the  nearest 
neighbors  of  the  defect  as  well,  are  reduced  as  the  overlap 
decreases.  This  allows  electronic  charge  to  flow  into  the  de¬ 
fect  region. 

For  interstitial  defects,  the  average  distance  between  the 
nearest  neighbors  and  the  atoms  further  out,  to  which  they 
are  bonded,  decreases  as  the  nearest  neighbors  relax 
outward.  This  increases  the  covalent  coupling  and  the 
charge  transfer  from  cations  to  anions  in  bonds  between  the 
nearest  neighbors  and  atoms  further  out.  In  order  to  mini¬ 
mize  fluctuations  due  to  this  charge  transfer  between  anions 
and  cations,  we  have  defined  an  interstitial  region  consisting 
of  the  interstitial  and  the  neighboring  atoms  which  are  close 
enough  interact  directly  with  the  interstitial  (i.e.,  the  first 
and  second  neighbor  shell  for  tetrahedral  interstitials,  and 
the  first  neighbor  shell  for  hexagonal  interstitials ) . ’  This  re¬ 
gion  contains  both  anion  and  cation  neighbors.  The  net 
amount  of  charge  transferred  into  the  interstitial  region 
when  the  lattice  is  allowed  to  relax  is  shown  in  Fig.  1  for  all 
tetrahedral  interstitials.  The  radial  relaxation  outward  is 
given  in  this  figure  as  the  fractional  change  in  the  interstitial- 
nearest  neighbor  distance  d .  A  clear  correlation  between  the 
amount  of  charge  transferred  and  the  preferred  amount  of 
relaxation  can  be  seen.  For  the  hexagonal  Cd  interstitial  in 
CdTe,  a  5.3%  increase  in  puckering  of  the  ring  of  nearest 
neighbors  is  allowed,  in  addition  to  the  preferred  6.3%  radial 
relaxation  of  the  nearest  neighbors  outward,'  resulting  in  a 
larger  charge  transfer  of  0.25  electrons  into  the  interstitial 
region.  For  the  cation  vacancy  in  HgTe,  there  is  also  a  small 
charge  transfer  (  -0.02  electrons)  onto  the  nearest  neigh¬ 
bors  of  the  vacancy  ,  resulting  from  their  2%  radial  relaxa¬ 
tion  inward. 

Large  overlaps  for  cation  intersiitials  between  cations 
(especially  Cd  between  Cd),  as  well  as  for  hexagonal  inter¬ 
stitials,  lead  to  larger  relaxations  outward  for  tetrahedral 
cation  interstitials  between  cations  than  for  cation  intersti¬ 
tials  in  the  preferred  tetrahedral  position  (between  anions). 
The  nearest  neighbors  of  the  Cd  interstitial  between  Cd  in 
CdTe  relax  outward  by  about  6%,  compared  with  relax¬ 
ations  of  3%-4%  for  Cd,  Hg,  and  In  interstitials  in  the  pre¬ 
ferred  tetrahedral  position  in  Hg,.,  CdrTe.' 
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Fig.  I.  Plot  of  predicted  minimum  energy  relaxation  versus  the  number  of 
electrons  transferred  into  the  interstitial  region  upon  relaxation. 


For  anion-like  interstitials  in  the  preferred  tetrahedral  po¬ 
sition,  most  of  the  electronic  charge  is  transferrecUo  the  in¬ 
terstitial  itself  as  the  lattice  relaxes.  During  relaxation,  less 
charge  is  transferred  to  the  anion-like  interstitials  which 
have  a  larger  excess  negative  charge  before  relaxation. 
Therefore,  we  can  predict  the  relative  amounts  of  relaxation 
outward  (3.1%,  4.6%,  and  5.2%  for  As.  Sb.and  Te)  simply 
by  knowing  the  amount  of  excess  negative  charge  on  these 
interstitials  before  relaxation:  0.26.  0.10,  and  0.01  e  for  As, 
Sb,  and  Te,  respectively,  in  units  of  the  electronic  charge. 

V.  CONCLUSIONS 

Results  of  self-consistent,  tight-binding  supercell  calcula¬ 
tions  for  the  localized  electronic  levels  and  lattice  relaxations 
associated  with  interstitial  As,  Sb.  and  Te  in  CdTe  and  the 
Hg  vacancy  in  HgTe  have  been  presented,  and  trends  m  the 
localized  levels  and  relaxations  for  these  defects  and  for  Cd, 
Hg,  and  In  interstitials,  from  earlier  calculations, 1  have  been 
discussed.  The  reduction  of  overlap  energies,  resulting  in  a 
net  transfer  of  electronic  charge  into  the  defect  region,  drives 
the  relaxation  for  all  these  defects.  Some  specific  conclusions 
which  may  be  drawn  from  these  results  are  given  below: 

Of  interstitials  in  the  preferred  tetrahedral  position,  Te 
has  the  largest  nearest-neighbor  relaxation  outward,  and  is 
the  most  likely  to-be  attracted  to  regions  where  this  extra 
strain  may  be  accommodated.  In  Te-rich  material,  there 
may  also  be  high  concentrations  of  cation  vacancies,  which 
have  an  inward  nearest-neighbor  relaxation  and  may  there¬ 
fore  attract  Te  interstitials.  This  strain-mediated  attraction 
may  promote  the  creation  of  Te  antisite  defects.  TheJarge 
strain  energy  associated  with  isolated  Te  interstitials  ma\ 
also  help  explain  the  tendency  for  excess  Te  to  condense  into 
second-phase  Te  inclusions. 

Large  strains  are  also  associated  with  Cd  interstitials  at 
hexagonal  sites  and  tetrahedral  sites  with  cation  nearest 


neighbors.  This  may  explain  the  apparent  unfavorability2  - 
of  Cd  interstitial  diffusion  down  the  open  channel  between 
preferred  tetrahedral  sites,  since  this  path  passes  through 
tetrahedral  sites  with  cation  nearest  neighbors  and  hexagon¬ 
al  sites. 

Defect  level  energies  for  anion-like  interstitials  in  their 
preferred  tetrahedral  position  (between  cations)  and  for  ca¬ 
tion-like  interstitials  in  their  preferred  tetrahedral  position 
(between  anions)  both  show  a  downward  trend  with  de¬ 
creasing  atomic  orbital  energies  on  the  interstitial. 

Vydyanath  has  suggested  that  the  midgap  tunneling 
centers  in  Hg-rich  HgCdTe  may  be  due  to  Hg  interstitials, 
based  on  his  experimental  work  in  various  overpressures  of 
Hg.18  Previous  calculations  using  the  less  accurate  small- 
basis  method  have  shown. a  deep  level  of  Ax  symmetry  in 
the  gap  region,  due  to  the  Hg  interstitial  in  the  preferred 
tetrahedral  position.  In  this  work,  we  find  that  the /I,  level 
which  holds  the  last  two  electrons  for  the  Hg  interstitial  in 
the  preferred  tetrahedral  position  appears  to  be  located  close 
to  the  conduction  band  edge  in  Hg-rich  HgCdTe.  This  inter¬ 
stitial  therefore  remains  a  candidate  fer  providing  tunneling 
levels  in  the  gap.  The  3%  outward  relaxation  of  the  Te  near¬ 
est  neighbors  of  this  interstitial  allows  th  e/I ,  level  to  become 
more  localized  on  the  interstitial  itself,  and  to  move  slightly 
lower  in  the  gap. 
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The  processes  of  formation  of  the  localized  defect  states  due  to  tetrahedral  and  hexagonal  sp~ type 
interstitials  in  sp-bonded  materials  are  considered  using  a  simple,  tight-binding-type  local-matrix  ap¬ 
proach.  In  this  method  we  subsequentK  account  for  the  interaction  of  the  interstitial  with  its 
nearest  neighbors,  and  for  the  covalent,  polar,  as  well  as  metallic  couplings  in  the  small  system  of 
bonds  surrounding  the  defect.  The  deep  and  hyperdeep  levels,  as  well  as  resonances  of  different 
symmetries,  are  discussed  and  are  shown  to  derive  either  from  the  atomic  states  on  the  interstitial, 
or  from  various  host  states  on  the  surrounding  atoms.  The  level-formation  schemes  for  the  P  inter¬ 
stitial  in  GaP,  self-interstitials  in  Si,  and  Hg  interstitials  in  CdTe  are  presented. 


I.  INTRODUCTION 


During  the  past  decade  considerable  progress  in  calcu¬ 
lations  of  localized  defect  states  has  been  made.  A 
variety  of  approaches  based  on  the  Green's-function 
method,  and  ranging  from  simple  semiempiricai  to  so¬ 
phisticated  self-consistent  local-density  formalisms,  has 
been  used  to  improve  our  knowledge  of  the  basic  proper¬ 
ties  of  localized  states.  However,  no  simple  interpreta¬ 
tion  of  the  results  obtained  follows  from  the  Green’s- 
function  method  itself.  In  particular,  neither  the  origins 
of  the  complicated  spectrum  of  localized  levels  and  reso¬ 
nances,  nor  the  character  of  these  levels,  can  be  explained 
by  the  method.  The  Green’s-function  approaches  also 
shed  no  light  on  the  important  role  of  the  local  environ¬ 
ment  of  the  defect  in  the  process  of  formation  of  localized 
states,  giving  no  transparent  picture  of  the  formation  of 
these  states  from  the  atomic  states  on  the  defect  and  the 
surrounding  host  atoms. 

To  interpret  the  Green’s-function  results,  simple,  usu¬ 
ally  defect-molecule-like,  local  models  have  been  used  for 
'substitutional1  and  interstitial2-4  defects.  However,  due 
to  their  often  schematic  character,  or  the  simplifying  as¬ 
sumptions  made,  the  presently  existing  interpretations 
give  a  limited  and  often  incomplete  physical  insight  into 
'  *iht  process  of  the  formation  of  localized  states. 

_  _  In  this  paper  we  give  a  far  more  detailed  analysis  of  the 
^formation  processes  of  localized  states,  focusing  on  sp - 
i^type  interstitials  in  tetrahedrally  bonded  semiconductors. 
p|We  use  the  local-matrix  approach  developed  in  the 
^preceding  paper5  (to  be  referred  to  as  I),  and  apply  it  to 
j^Etctrahedral  and  hexagonal  interstitials.  As  discussed  in  I, 
k^nis  simple  tight-oinding-type  approach  results  in  a  good 
Qualitative  and  often  quantitative  agreement  with  more 
jiophisticated  self-consistent  Green’s-function  calcula- 
jon s,  and  thus  may  provide  us  with  a  fairly  realistic  pic- 
je  of  the  formation  processes.  The  method  allows  us  to 
ider  the  formation  of  the  localized  states  not  only  in 
opolar  crystals  such  as  Si,2-4  but  also  in  materials 
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with  strong  polar  and  metallic  couplings  such  as 
Hg^Cd.Te.  We  interpret  both  deep  and  hyperdeep 
states,  as  well  as  hosthke  resonances  (or  levels),  often 
present  in  the  midregion  of  the  valence  band  of 
tetrahedrally  bonded  semiconductors,  whose  origins  have 
not  yet  been  analyzed.  The  distributions  of  charge  densi¬ 
ty  associated  with  the  localized  states  are  also  directly  in¬ 
terpreted  in  terms  of  the  parent  states. 

In  Sec.  II  we  briefly  present  the  local-matrix  method. 
Then,  following  closely  the  course  of  discussion  adopted 
in  paper  I,  in  Sec.  Ill  we  examine  the  formation  processes 
for  the  tetrahedral  P  interstitial  in  GaP,  for  tetrahedral 
and  hexagonal  self-interstitials  in  Si,  and  finally  for  Hg 
tetrahedral  and  hexagonal  interstitials  in  CdTe. 


II.  METHOD 

Let  us  briefly  summarize  the  basic  assumptions  of  the 
local-matrix  method,5  which  sets  a  framework  for  the 
present  analysis  of  formation  mechanisms.  In  this 
method,  a  local  system  of  bonds  surrounding  the  defect  is 
softly  decoupled  from  the  infinite  jp-bonded  host  crystal, 
and  described  in  terms  of  the  local  matrix.  The  soft- 
decoupling  procedure  preserves  all  the  covalent  (strong¬ 
est)  bonds,  while  neglecting  other,  weaker  interactions 
that  couple  the  local  system  to  the  rest  of  the  host  crys¬ 
tal.  In  addition  to  the  covalent  interaction,  we  account 
for  the  metallic  coupling  (promotion  energy)  and  polar 
coupling  in  the  system.  The  interaction  between  the  in¬ 
terstitial  and  its  nearest  neighbors  is  treated  exactly, 
within  the  assumed  two-center  approximation.6  The  in¬ 
teraction  between  the  host  atoms  is  considered  within  the 
nearest-neighbor  approximation.  The  values  of  tight- 
binding  parameters  entering  the  local  matrix,  as  well  as 
the  reference  position  of  the  valence-band  edge,  are  the 
same  as  in  the  host  crystal  from  which  the  local  system 
has  been  decoupled.  The  position  of  the  top  of  the 
valence  band  is  thus  given  by  the  sp 3  nearest-neighbor 
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tight-binding  calculation  for  the  perfect  host  crystal,  with 
the  same  values  of  tight  binding  parameters  as  those  used 
m  the  local  matrix. 

For  the  tetrahedral  and  hexagonal  interstitials  we  con¬ 
sider  the  locc*l  matrix  that  involves  36  and  40  basis  states 
(s  and  p  orbitals,  and/or  sp}  hybrids)  located  on  15  and 
19  atoms,  respectively.  Namely,  we  include  in  the  calcu¬ 
lation  all  states  on  the  interstitial  and  on  its  nearest 
neighbors  (INN  atoms),  and  all  those  hybrids  on  further 
neighbors,  which  are  covalently  bonded  to  hybrids  on  the 
nearest  neighbors  (the  further-neighbor  atoms  that  have 
covalent  bonds  with  the  INN  atoms  will  be  called  the 
2NN  atoms;  see  Fig.  1  in  I).  The  values  of  the  on-site 
tight-binding  matrix  elements  E'*}  and  Ep*\  where  k 
denotes  the  atomic  species,  are  taken*  from  atomic  data 
for  the  s  and  p  levels,  respectively.  The  off-site  tight- 
binding  matrix  elements  for  both  host-host  and  host- 
defect  nearest-neighbor  interactions  assume  the  empirical 
1  /d2- scaled  values.7 

At  each  stage  of  the  formation  process  we  solve8  the  ei¬ 
genvalue  equation  for  the  local  matrix  considered  in  I,  in 
which  appropriate  (for  the  considered  stage)  tight- 
binding  interactions  are  taken  into  account.  We  start 
with  the  system  composed  of  isolated  atoms  (i.e.,  the  sys¬ 
tem  with  all  off-site  elements  neglected),  and  subsequently 
include  the  interactions  between  the  atoms  in  order  of  the 
increasing  importance.  To  characterize  the  states  of  the 
system  at  each  step  of  the  formation  process,  along  with 
the  energy  of  each  state,  we  also  calculate  the  corre¬ 
sponding  probability  of  finding  the  electron  on  the  inter¬ 
stitial  site,  on  all  INN  atoms,  and  on  all  2NN  atoms. 
This  probability  (or  charge)  distribution,  expressed  in 
percent,  is  given  in  parentheses  for  each  energy  level 
shown  in  the  figures  discussed  in  the  next  section.  All  pa¬ 
rameters  used  in  the  present  calculation  assume  the 
values  listed  in  Sec.  V  of  I. 


III.  RESULTS 

In  this  section  we  discuss  the  formation  mechanisms 
for  several  tetrahedral  and  hexagonal  interstitials  in  GaP, 
Si,  and  CdTe,  following  the  course  of  discussion  present¬ 
ed  in  I. 


A.  The  tetrahedral  P  interstitial  in  GaP 

The  localized  states  obtained  in  the  local-matrix  ap¬ 
proach  for  the  neutral  P  interstitial  (with  P  nearest  neigh¬ 
bors)  at  the  tetrahedral  (Td  symmetry)  site  in  GaP  have 
been  discussed  in  I,  giving  a  good  agreement  with  the  re¬ 
sults  of  self-consistent  Green’s-function  calculations.9 
Here  we  will  consider  the  process  of  formation  of  these 
states,  winch  is  presented  in  Fig.  1,  and  consists  of  the 
following  stages. 

(i)  Panel  (a)  shows  the  energy  levels  E{sb)  and  E[pb)  of  the 
isolated  interstitial  (k=8)  atom,  and  the  energy  level  E 5 
of  the  dangling  sp 3  hybrids  on  the  INN  atoms.  Here, 
Ek  =z~(E{JK)jr3EpK))>  with  *=1,2  for  the  INN  and  2NN 
atoms,  respectively.  All  the  interactions  between  orbitals 
on  different  atoms,  and  the  promotion  energy  on  the 
INN  atoms,  are  neglected  at  this  stage.  The  states  of  the 


nomntei  acting  system  corresponding  to  energies  E,  )  and  , 
L'p'  have  all  the  charge  located  on  tne  interstitial,  and 
the  states  with  the  energy  Ex  have  all  the  charge  distri¬ 
buted  among  the  INN  atoms.  We  observe  that  the  sys¬ 
tem  composed  of  noninteracting  atoms  also  has  states 
with  energy  E 2,  t.e.,  the  energy  of  atomic  states  on  the 
2NN  atoms  twe  note  that  according  to  the  soft- 
decoupling  procedure  we  always  neglect*  the  promotion 
energy  on  the  2NN  atoms).  These  states,  whose  charge  is 
distributed  entirely  among  the  2NN  atoms,  do  not  parti¬ 
cipate  in  the  first  stage  of  formation,  and  will  be  included 
later. 

Ui)  Panel  (b)  shows  the  splitting  of  the  levels,  which  is 
due  to  the  interaction  of  the  dangling  hybrids  on  each 
INN  atom  with  the  interstitial  hybrid  that  points  toward 
this  atom.  [This  interaction  is  described10  by  the  matrix 
elements  ps±  and  ppT  given  by  Eqs.  (13)  and  (14)  of  I,  re- 
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FIG.  1.  Formation  of  the  localized  states  due  to  the 
tetrahedral  P  interstitial  (with  P  nearest  neighbors)  in  GaP. 
Formation  stages:  (a)  isolated  atoms;  (b)  including  the  (leading 
part  of  the)  interaction  of  the  interstitial  orbitals  with  dangling 
hybrids  on  the  INN  atoms;  (c)  including  the  covalent  coupling 
between  hybrids  on  the  INN  and  2NN  atoms;  (d)  including  the 
promotion  energy  on  the  INN  atoms;  (e)  including  the  remain' 
ing  interactions  between  the  interstitial  and  the  INN  atoms  (see 
details  in  text).  Notation:  localized  states,  thick  solid  lines* 
weakly  localized  or  extended  states,  dashed  lines;  parent  host 
states,  thin  solid  lines;  la,6,c),  fractional  charge  (in  percent)  c 
the  interstitial,  on  all  INN  atoms,  and  on  all  2NN  atoms,  re-  ^ 
spectively.  The  energies  of  the  system  without  the  intcrstiti 
EJ±  and  Ep±  (see  text)  are  indicated  on  the  right-hand  side  < 
the  figure. 
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■  spectivelv.]  The  rest  of  the  defect-host  interactions,  less 
significant  for  the  formation  process,  will  be  included 
later.  The  levels  split  into  a  pair  of  bondinglike  levels  of 
A  i  (s-hke)  and  7\  [pA ike)  symmetry  located  in  the 
valence  band,  and  a  pair  of  A  j  and  7\  antibondinghke 
levels  located  above  the  top  of  the  valence  band. 

dii)  Panel  (c)  shows  the  result  of  including  the  covalent 
coupling  /?  between  the  host  atoms;  1  e  ,  the  dangling  h>- 
brids  on  the  INN  atoms  now  form  covalent  bonds  with 
the  hybrids  on  the  2NN  atoms.  Due  to  this  interaction, 
the  three  lower  levels  of  panel  (b)  are  shifted  downward  in 
energy,  while  the  upper  T2  level  goes  upward,  loses  a  lot 
of  charge  on  the  interstitial  site,  and  becomes  a  weakly 
localized  level  (denoted  by  a  dashed  line)  A  pair  of  new 
A  |  and  T 2  states  is  also  formed.  The  formation  of  these 
new  states  may  be  seen  as  resulting  from  the  interaction 
of  the  already  existing  states  tof  A  i  and  T2  symmetry,  re¬ 
spectively)  with  the  above-mentioned  host  states  (of  the 
system  composed  of  noninteracting  atoms)  localized  on 
the  2NN  atoms  and  having  the  energy  £\.  The  energy 
level  corresponding  to  these  states  is  indicated  in  panel 
(b)  at  about  4*4  eV.  Due  to  this  interaction  all  the  previ¬ 
ously  existing  levels  are  shifted  as  described  above,  while 
the  E2  level  splits  into  a  pair  of  a  strongly  localized  7% 
level  and  a  delocalized  A ,  level  (denoted  by  a  dashed 
line).  We  observe,  however,  that  because  the  E:  level  and 
the  T2  level  at  +4  eV  in  panel  (b)  are  very  close  to  each 
other,  the  resulting  localized  7\  state  [at  about  4-2  eV  in 
panel  (c)]  has  a  strongly  mixed  character.  Hence,  its  as¬ 
sociation  with  the  parent  level  £2  has,  m  this  case,  a  rath¬ 
er  conventional  character,  and  we  may  as  well  consider 
this  T2  level  to  be  derived  from  the  T2  level  at  4-4  eV  in 
panel  (b). 

i‘  (iv)  Panel  (d)  shows  the  result  of  including  the  promo¬ 
tion  energy  At  on  the  INN  atoms,  i.e.,  now  the  hybrids 
'on  the  same  INN  atoms  can  interact  with  each  other; 
^here,  A,  =  —  ±(E}n—Epl>).  For  the  states  of  Al  symme¬ 
try,  the  most  significant  effect  due  to  this  interaction  is 
“the  formation  of  a  third  localized  state  of  A  {  symmetry, 
located  below  the  top  of  the  valence  band.  The  formation 
bfthe  third  localized  A  t  level  may  be  viewed  as  resulting 
from  the  interaction  of  the  two11  already  existing  A{  lo¬ 
calized  states  with  the  bonding  A  {  symmetry  host  state 
of  the  system  formed  due  to  the  covalent  (and  also  polar) 
Interaction  between  the  INN  and  2NN  atoms  (namely, 
he  A  {  state  of  the  system  without  a  defect,  with  the  en- 
*gy  given  by  Eq.  (5)  of  I,  taken  for  the  minus  sign,  and 
:h  A^O,  i.e.,  with  the  energy  E^=y(£l 
[(£,  —  £2)24-4/?2]1/2  j ).  This  host  state,  with  the 
?rgy  £_,  is  located  between  the  two  A{  localized 
Wes  of  panel  (c),  close  to  —4  eV.  The  interaction  be- 
geen  ihese  three  states  pushes  the  upper  A ,  localized 
fel  towards  the  top  of  the  valence  band  and  pushes 
the  lower  A  {  level  out  of  the  valence  band,  while 
host  state  is  pulled  down  into  the  valence  band,  as- 
the  localized  character,  and  forms  the  third  A ,  Io- 
level.  A  similar  process  occurs  for  the  T2  states. 
,  we  have  two  T2  localized  states  interacting  with 
bonding  T2  symmetry  host  states  formed  due  to  the 
interaction  between  the  INN  and  2NN  atoms 
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[i.e.,  the  T2  states  of  the  system  without  a  defect,  with  the 
energy  given  by  Eq.  lb)  of  I,  taken  for  the  minus  sign,  and 
with  A,  -0].  The  host  level  (with  the  same  energy  £-  as 
in  the  case  of  the  A ,  host  state,  considered  above)  is 
again  located  between  the  two  localized  T2  levels  of  panel 
lcj.  Due  to  the  interaction  of  these  three  states,  the  upper 
7\  localized  level  remains  almost  unaffected,  the  host  lev¬ 
el  is  slightly  pushed  down  in  energy,  assumes  a  localized 
character,  and  forms  the  second  T2  level  [at  around  —4 
eV  in  panel  la)],  while  the  lower  7\  ^evei  strongly 
pushed  down  in  energy,  loses  a  lot  of  charge  on  the  inter¬ 
stitial  site,  and  becomes  a  weakly  localized  level  (denoted 
by  a  dashed  line).  We  note  that  inclusion  of  the  promo¬ 
tion  energy  on  the  INN  atoms  leads  also  to  the  interac¬ 
tion  with  the  antibonding  Ax  and  T2  symmetry  host 
states  of  the  system  with  the  energy  £.*.,  formed  due  to 
the  covalent  interaction  between  the  INN  and  2NN 
atoms  (namely,  the  A  j  and  T2  states  of  the  system 
without  a  defect,  with  energies  given  by  Eqs.  (5)  and  (6)  of 
I,  respectively,  taken  for  the  plus  sign,  and  with  A|=0, 
i.e.,  with  the  energy  £  +  =  t{£i  +E2  4-[(£j 
~£2  )24-4/?2]i/2  j ),  and  located  at  around  4-6  eV.  How¬ 
ever,  this  coupling  results  in  states  that  are  not  localized, 
and  thus,  for  the  sake  of  clarity,  not  displayed  in  Fig.  1. 

tv)  Panel  (e)  shows  the  result  of  including  the  remain¬ 
ing  interactions  between  the  interstitial  and  the  INN 
atoms  [i.e.,  those  described  by  the  matrix  elements  v±, 
7?±,  and  p±  given  by  Eqs.  (15M17)  of  1].  The  main  effect 
due  to  these  interactions  is  the  further  reduction  of  the 
Ax  —  T2  splitting  of  the  levels  in  the  fundamental  gap  re¬ 
gion,  and  further  lowering  of  the  hyperdeep  A  \  level.  On 
the  right-hand  side  of  Fig.  1  we  also  show,  for  compar¬ 
ison,  the  energies  Es±  and  Epx  [given  by  Eqs.  (5)  and  (6) 
of  I]  of  the  system  without  the  interstitial  (but  with  all 
other  interactions  included).  To  each  of  these  levels  there 
corresponds  one  state  of  A  x  symmetry  and  three  states  of 
7\  symmetry. 

The  above  analysis  indicates  that  two  of  the  obtained 
localized  states  may  be  seen  as  originating  from  the  atom¬ 
ic  states  of  the  interstitial;  namely,  the  T2  deep  state  in 
the  gap  derives  from  the  p  atomic  state  of  the  interstitial 
isee  the  discussion  above)  and  the  lowest  A  i  hyperdeep 
states  derives  from  the  s  atomic  state  of  the  interstitial. 
The  large  amount  of  charge  on  the  interstitial  found  for 
these  states,  as  well  as  a  strong  variation  of  their  energies 
with  changes  in  the  atomic  energies  on  the  interstitial, 
art  consistent  with  this  picture.  Therefore,  these  two 
states  may  be  called  defectlike.  We  note  that  both  of 
these  states  are  substantially  mixed  with  the  host  states; 
for  example,  the  deep  Tz  state  is  strongly  mixed  with  the 
Ez  state  localized  on  the  2NN  atoms,  and  because  of  that 
it  has  most  ot  the  remaining  charge  localized  on  the  2NN 
atoms. 

The  remaining  three  localized  states  have  a  hostlike 
character:  the  gap-region  A ,  localized  state  derives  from 
the  host  states  of  the  INN  atoms,  while  a  pair  of  A  x  and 
T2  mid-valence-band  states  derive^  from  the  host  states 
formed  from  covalent  bonds  between  the  INN  and  2NN 
atoms.  All  of  these  three  states  have  smaller  charges  on 
the  interstitial  and  weaker  variations  of  their  energies 
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wnh  the  change  of  atomic  energies  on  the  interstitial. 
These  variations  of  energies  are  especially  weak  for  the 
mid-valence-band  states,  which  reflects  their  strongly 
hostlike  character,  and  are  consistent  with  the  behavior 
observed  for  such  states  m  the  self-consistent  Green's- 
function  calculations. 

B.  Tetrahedral  self-interstitials  in  Si 

We  shall  now  consider  the  formation  of  localized  levels 
due  to  the  neutral  tetrahedral  self-interstitial  in  Si  The 
comparison  of  the  results  of  the  present  approach  with 
the  self-consistent  Green’s-function  calculations12  have 
been  discussed  in  I.  Since  many  features  of  the  formation 
process,  presented  in  Fig.  2,  are  similar  to  the  previously 
considered  case  of  the  P  interstitial  in  GaP,  in  our  discus¬ 
sion  below  we  will  mainly  focus  on  differences  between 
these  two  cases. 

As  for  the  P  interstitial  in  GaP,  the  £,  level  corre¬ 
sponding  to  host  states  located  on  the  INN  atoms  lies  be¬ 
tween  the  E'b)  and  Epb)  energies  on  the  interstitial,  and 
the  levels  resulting  due  to  the  leading  part  of  the  interac¬ 
tion  between  the  interstitial  and  the  INN  atoms,  shown 
in  panel  lb),  are  qualitatively  similar  to  those  obtained  for 
P  in  GaP.  The  difference  occurs  when  the  covalent  cou¬ 
pling  ft  between  the  INN  and  2NN  atoms  is  included. 
Now  the  energy  £2  (^2  =  £ , )  of  the  host  states  located 
on  the  2NN  atoms  lies  between  the  two  pairs  of  A  {  and 
Tl  levels  of  panel  (b).  Due  to  this,  both  the  and  f2 
states  formed  in  the  gap  region  [panel  (c)]  derive  form  the 
host  £2  states,  and  have  most  of  the  charge  located  on 
the  2NN  atoms.  Now,  unlike  for  P  in  GaP,  also  the  hy- 


FIG.  2.  Formation  of  the  localized  states  due  to  the 
tetrahedral  self-interstitial  in  Si.  Formation  stages  and  notation 
as  in  Fig.  1. 
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perhigh  7\  level  (resonance)  in  the  conduction  band  is  lo¬ 
calized.  The  inclusion  of  the  promotion  energy  on  the 
INN  atoms  [panel  (d)]  results,  similarly  as  for  P  in  Gap, 
in  the  formation  of  a  pair  of  new  A  j  and  7\  states  locat¬ 
ed  in  the  middle  of  the  valence  band.  They  derive  from 
the  £_  (E±  =£|±i/3|)  bonding  Ax  and  7%  symmetry 
host  states  of  the  system,  respectively  (see  Sec.  Ill  A  for 
details/.  We  note  that  the  £„  states  now  contain  more 
charge  on  the  2NN  atoms  than  in  the  case  of  P  in  GaP, 
and  due  to  that  the  localization  of  the  resulting  A  {  state 
tat  about  —  5  eV)  is  weaker  than  for  the  comparable  state 
for  P  in  GaP,  while  the  localization  of  the  T2  state  (at 
about  —4.5  eV)  becomes  very  weak,  and  thus  this  level  is 
denoted  by  a  dashed  line.  As  in  the  case  of  the  P  intersti¬ 
tial  in  GaP,  the  antibonding  host  states  with  the  energy 
E+  do  not  give  rise  to  any  localized  states.  The  inclusion 
of  the  remaining  direct  interactions  between  the  intersti¬ 
tial  and  the  INN  atoms  [panel  (e)j  leads  to  the  same 
effects  as  for  P  in  GaP.  We  observe  that,  similar  to  that 
case,  these  interactions  do  not  practically  affect  the  posi¬ 
tions  of  the  levels  of  7%  symmetry.  The  distribution  of 
the  total  charge  between  the  interstitial,  all  the  INN 
atoms,  and  all  the  2NN  atoms  is  also  weakly  affected  by 
these  interactions,  though  the  wave  function  itself  can  ex¬ 
perience  significant  changes  (up  to  100%)  on  single  sites. 
It  is  also  worthwhile  to  note  that  although  the  A  x  state 
located  at  about  —  2  eV  in  panel  (e)  has  most  of  its  charge 
located  on  the  2NN  atoms,  it  is  still  a  well-localized  state. 
It  has  29%  of  charge  located  on  the  interstitial,  about 
3%  on  each  INN  atom,  about  7%  on  each  2NN  atom 
from  the  second  shell,  and  about  5%  on  each  2NN  atom 
from  the  fourth  shell.  Thus,  the  charge  on  the  interstitial 
site  is  significantly  larger  than  on  any  other  single  site, 
which,  according  to  our  criteria  discussed  in  I,  classifies 
this  state  as  localized. 

The  above  analysis  of  the  level  formation  for  the 
tetrahedral  self-interstitial  in  Si  shows  that,  as  in  the  case 
of  P  in  GaP,  the  hyperdeep  A  \  state  has  a  defectlikc 
character,  deriving  from  the  s  atomic  state  on  the  inter¬ 
stitial.  We  note  the  strong  mixing  of  this  state  with  the 
£,  and  £_  host  states.  The  T2  gap  state  derives  from 
the  host  states  localized  on  the  2NN  atoms.  It  has,  how¬ 
ever,  a  very  strong  component  coming  from  the  Ep  state 
on  the  interstitial.  This  T2  state  also  has  a  large  charge 
on  the  interstitial,  and  the  strongest  variation  of  its  ener¬ 
gy  with  changes  in  the  E{pb)  energy  on  the  interstitial* 
Thus,  we  would  characterize  this  level  as  predominantly 
defectlike.  The  two  upper  A ,  localized  states  in  the 
valence  band,  as  well  as  a  weakly  localized  T2  state  at  ^ 
about  —4.5  eV,  are  hostlike.  They  have  smaller  charges 
on  the  interstitial,  and  weaker  variations  of  the  energies 
with  changes  of  atomic  energies  on  the  interstitial. 


C.  The  hexagonal  self-interstitial  in  Si 


The  process  of  formation  of  the  localized  states  due  t  • 
the  neutral  hexagonal  self-interstitial  in  Si  is  shown  . 
Fig.  3.  The  formation  stages  are  similar,  but  not  identi  ^ 
to  those  discussed  above  for  the  tetrahedral  interstitial*  £ 
and  are  as  follows. 
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FIG  3  Formation  of  the  localized  states  due  to  the  hexago¬ 
nal  self-interstitial  in  Si.  Formation  stages;  (a)  isolated  atoms; 
(b)  including  the  interaction  of  the  interstitial  orbitals  with 
dangling  hybrids  on  the  INN  atoms;  (c)  including  the  covalent 
coupling  between  hybrids  on  the  INN  atoms  in  the  hexagonal 
ring;  Id)  including  the  covalent  coupling  between  hybrids  on  the 
INN  and  2NN  atoms;  (e)  including  the  promotion  energy  on  the 
INN  atoms.  Notation  as  in  Fig.  1. 
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rstitial 


(i)  Panel  (a)  show  the  energy  levels  E(:b)  and  E{p6)  of  the 
isolated  interstitial  atom,  and  the  energy  level  Eh  of  the 
dangling  sp}  hybnds  on  the  INN  atoms.  Note  that  in  the 

^homopoiar  material  the  hybrid  energies  on  all  INN 
l  atoms  are  the  same:  Eh—E{=E2 .  All  interactions  be- 
tween  orbitals  on  different  atoms,  and  the  promotion  en¬ 
ergy  on  the  INN  atoms,  are  neglected  at  this  stage. 

(ii)  Panel  (b)  shows  the  splitting  of  the  levels  due  to  the 
,  ^eracti°ns  of  s  and  p  orbitals  on  the  interstitial  with  the 
[^dangling  hybrids  on  the  INN  atoms.  Due  to  the  lower 

symmetry  of  the  hexagonal  site,  the  meaningful 
['^fkcomposition  of  the  interstitial-host-atoms  interaction 
ato  more  and  less  significant  parts,  utilized  in  the 
^  ctrahedral  case,  cannot  be  made.  Therefore,  the  full  in- 
rg&teraction  between  the  interstitial  and  the  INN  atoms  is 
lincluded  here.  The  El  level  snlits  into  n  spt  of  honHina- 
Is  of  A  *  (s-like),  £“  (pXJpy- like),  and  A2  (p2- 
Ice)  symmetry  located  below  the  top  of  the  valence  band, 
°d  a  similar  set  of  antibondinghke  levels  located  in  the 
aduction  band.  The  px  and  py  orbitals  on  the  intersti- 
lying  in  the  plare  of  the  puckered  hexagon,  interact 
the  INN  atoms  more  strongly  than  the  pz  orbitals, 
1  leads  to  a  larger  splitting  between  the  bonding 
^antibonding  states  than  seen  for  the  A  7  states. 
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(iii)  Panel  (c)  shows  the  effect  of  including  the  covalent 
coupling  between  the  INN  atoms  that  form  the  hexago¬ 
nal  ring  around  the  interstitial:  now  two  hybnds  on  each 
INN  atom  are  involved  in  covalent  bonds  in  the  ring, 
while  two  remaining  hybrids  on  each  INN  atom  are  still 
dangling  (compare  to  Fig.  1  in  I).  The  coupling  strongly 
shifts  down  the  energies  of  A  {  and  hyperdeep  as 
well  as  hyperhigh  states.  Some  new  localized  states  are 
also  formed.  They  arise  from  the  Eh  host  states  located 
on  the  INN  atoms.  (The  degeneracy  of  the  Eh  level  was 
not  completely  removed  by  the  host-defect  interactions 
considered  in  the  previous  stage,  and  now  it  splits  again 
due  to  the  covalent  interactions  in  the  ring.13)  The  Eh 
level  splits  into  A  7  and  £“  localized  levels  placed  in  the 
conduction  band,  and  a  delocalized  A  {  state  (mostly 
composed  of  the  dangling  hybrids  on  the  INN  atoms) 
placed  just  below  the  top  of  the  valence  band.  We  ob¬ 
serve  a  strorg  shift  in  the  characters  of  the  two  A  7  levels 
placed  in  the  conduction  band.  The  upper  one  (at  +5 
eV)  is  strongly  localized  on  the  INN  atoms,  and  the 
lower  one  (at  4-3.5  eV)  has  a  significant  amount  of  charge 
on  the  interstitial.  However,  the  characters  of  the  levels 
they  derive  from  are  reversed,  namely  the  upper  A  7 
state  in  panel  (b)  is  strongly  localized  on  the  interstitial, 
while  the  lower  one  is  the  Eh  host  state  localized  on  the 
INN  atoms.  Since  the  shift  in  the  character  of  the  states 
is  so  drastic  here,  we  have  chosen  to  connect  (by  dashed 
lines)  the  A  2  levels  of  an  alike  character,  rather  than  to 
follow  the  movement  of  the  levels  when  (3  infractions  in 
the  ring  were  slowly  increased  from  zero  (in  the  latter 
case,  clearly,  the  levels  would  not  cross).  In  f  cuture 
cases,  when  there  is  such  a  clear  shift  in  the  character  of 
the  states  that  we  must  show  the  two  states  of  the  same 
symmetry  crossing  in  order  to  follow  the  character  of  the 
states,  we  will  use  similar  dashed  lines. 

(iv)  Panel  (d)  shows  the  effect  of  including  the  remain¬ 
ing  covalent  interactions  between  the  host  atoms,  i.e.,  the 
covalent  coupling  between  the  INN  and  2NN  atoms. 
[Note  that  panel  (d)  for  the  hexagonal  case  is  equivalent 
to  panel  (c)  for  the  tetrahedral  case.)  The  hyperdeep  A  f 
and  £“  states  are  almost  unaffected  by  these  interac¬ 
tions,  since  in  stage  (c)  these  states  have  almost  all  of  the 
charge  located  on  the  interstitial  or  engaged  in  the  ring 
bonds.  In  contrast,  the  lowest  A2  state  at  —2  eV  in 
panel  (c)  consists  mostly  of  a  bondinglike  combination  of 
dangling  hybrids  on  the  INN  atoms  and  tli% .  z  orbital  on 
the  interstitial.  Since  this  coupling  is  cons’derably  weak¬ 
er  than  the  covalent  bonds  between  die  INN  and  2NN 
atoms,  the  considered  A  2  state  loses  a  lot  of  charge  on 
the  interstitial  when  these  bonds  are  formed,  urtd  be¬ 
comes  delocalized.  This  charge  cn  the  interstitial  site 
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the  bottom  of  the  conduction  band  in  panel  (d).  For  the 
A  ,+  levels  in  the  gap  region  we  find  again  a  strong  shift  in 
the  characters,  which  is  why  these  levels  were  connected 
by  the  dashed  lines. 

(v)  Panel  (e)  shows  the  effect  of  including  the  promo¬ 
tion  energy  on  the  INN  atoms.  As  m  the  case  of  the 
tetrahedral  self-interstitial,  this  interaction  lowers  the  en¬ 
ergy  of  the  hyperdeep  states  and  increases  the  energy  of 
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the  hyperhigh  states.  Also,  similar  to  the  tetrahedral 
case,  it  leads  to  the  formation  of  a  new  A  j*  localized 
hostlike  state,  which  derives  from  the  bonding  E_ 
vEt  =  E,  ±,/3| )  host  states,  formed  due  to  the  covalent  in¬ 
teraction  between  the  INN  and  2NN  atoms  in  the  hexag¬ 
onal  system  with  no  interstitial. 

When  the  promotion  energy  on  the  INN  atoms  is  in¬ 
cluded,  some  states  derive  also  from  the  antibonding  host 
states  with  the  energy  E  ^  located  at  about  -*-4.5  eV>.  In 
the  case  of  the  tetrahedral  self-interstitial  the  resulting 
state  is  extended,  while  the  77  state  has  a  clear  sur¬ 
face  character.  However,  for  the  hexagonal  self- 
interstitial,  the  E  +  host  states  give  rise  to  the  A  {  state 
of  a  rather  weak  surface  character,  but  with  a  consider¬ 
able  amount  of  charge  labout  20%}  on  the  interstitial, 
and  the  energy  strongly  shifted  downward  (to  about  -r-  2.5 
cV).  This  could  suggest  that  we  have  obtained  a  second 
localized  A  7  state  in  the  gap  region.  However,  obvious¬ 
ly,  the  considered  A  7  hexagonal  state  splits  (when  going 
from  the  Td  to  the  D}d  site)  from  the  above-mentioned 
T2  surface  State.  This  strongly  indicates  that  the  con¬ 
sidered  A  2  state  is  a  spurious  quasilocalized  state  whose 
properties  derive  from  the  small  size  and  particular  form 
of  the  hexagonal  system,  i.e.,  such  a  state  should  not  be 
observed  as  localized  m  the  bulk  crystal.  This  conclusion 
is  supported  by  other  calculations,14  which  also  find  a  sin¬ 
gle  A  2  state  in  the  gap  region.  The  spurious  A  7  state 
formed  just  above  the  localized  A  7  level  may  be  seen  as 
being  responsible  for  the  incortect  ordeimg  of  levels  seen 
for  the  hexagonal  Si  self-mterstitiai  in  the  gap  region:  the 
A  f  level  should  lie  below1415  the  A  7  and  E“  levels,  as 
it  does  for  the  other  three  groups  of  states  in  panel  (e). 
Namely,  the  spurious  A  2  state,  which  acquires  a  lot  of 
charge  on  the  defect  site  from  the  localized  A  7  gap  state, 
repels  it  further  down  in  energy,  which  is  why  the  A  2  lo¬ 
calized  state  stays  below16  the  A  f  gap  state. 

Comparing  the  final  structures  of  the  levels  obtained 
for  the  tetrahedral  and  the  hexagonal  self-interstitials,  we 
find  that  (as  it  should  be)  to  each  tetrahedral  A ,  level 
there  corresponds  one  hexagonal  A  ?  level,  while  each 
tetrahedral  T2  level  splits  into  a  pair  of  A  7  and  £”  hex¬ 
agonal  levels,  with  the  characters  of  the  corresponding 
levels  basically  preserved  during  the  transition  from  the 
Td  to  the  D}d  site.  We  observe  that  when  this  transition 
occurs,  the  very  weakly  localized  T2  hyperdeep  level 
gives  rise  to  the  localized  E_  hyperdeep  level,  seen  also 
in  other  calculations.14,15 


D.  Tetrahedral  and  hexagonal  Hg  interstitials  in  CdTe 

In  paper  I,  we  calculated  the  localized  states  due  to  in¬ 
dium  and  self-interstitials  in  Hg^Cd^Te  Here,  as  an 
example,  we  have  chosen  to  discuss  the  formation  of  the 
localized  states  due  to  the  Hg  interstitial  in  pute  CdTe. 
In  the  tetrahedral  case,  we  will  consider  only  the  pre¬ 
ferred  position  of  *he  Hg  interstitial,  i.e.,  that  with  anion 
nearest  neighbors;17  this  interstitial  may  be  responsible 
for  providing  the  midgap  tunneling  levels  in  the 
mercury-nch  compound.5  We  note  that  since  In,  and 
also  Cd,  have  values  of  the  s  and  p  atomic  energies  close 
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to  those  of  Hg,  the  results  shown  belo'*  are  also  typical 
for  these  interstitials. 

Let  us  start  with  the  formation  process  due  to  the 
tetrahedral  Hg  interstitial  (with  Te  nearest  neighbor:' 
shown  in  Fig.  4.  We  first  observe  that,  unlike  in  the  pre¬ 
viously  considered  case  of  P  in  GaP  (i.e.,  also  an  intersti¬ 
tial  with  anion  nearest  neighbors),  or  in  the  case  of  the 
self-interstitial  in  Si,  now  both  and  E{pb)  energies  on 
the  defect  lie  above  the  E,  energy  on  the  INN  atoms 
[panel  (a)].  This  change  of  the  ordering  of  the  parent  lev¬ 
els  on  isolated  atoms  has  profound  implications  on  the 
character  of  some  finally  formed  localized  states.  In  par¬ 
ticular,  the  finally  formed  A ,  gap  state  [at  about  4- 1  eV 
in  panel  (e)]  now  has  a  defectlike  character,  deriving 
directly  from  the  atomic  s  state  on  the  interstitial  and 
having  half  of  its  charge  located  on  the  interst/ial. 
Furthermore,  the  hyperdeep  A  l  state  is  now  delocalized, 
while  in  the  previously  considered  cases  the  comparable 
A ,  hyperdeep  states  had  a  strongly  localized  defective 
character.  The  mid-valence-band  A ,  resonance  is,  as  in 
the  previously  considered  formation  processes,  of 
predominantly  hosthke  character,  with  most  of  its  charge 
located  on  the  INN  atoms.  Among  the  states  of  77  sym¬ 
metry  the  lowest  localized  one  is  placed  in  the  conduction 
band  at  about  4  4  eV.  This  predominantly  hostlike  state 
originates  from  atomic  states  on  the  2NN  atoms,  and  has 
most  of  its  charge  located  on  the  2NN  atoms  from  the 
second  shell,  with  only  minute  charges  on  the  INN  atoms 
and  on  the  2NN  atoms  from  the  fourth  shell.  We  may 
observe  a  similarity  of  the  character  of  this  77  state  to 
the  77  gap  states  found  previously  for  P  in  GaP  and  for 
the  self-interstitial  in  Si.  The  hyperhigh  77  state  now  has 


FIG.  4.  Formation  of  the  localized  states  due  to 
tetrahedral  Hg  interstitial  (with  Te  nearest  neighbors)  in 
Formation  stages  and  notation  as  in  Fig  i 
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a  defectlike  character,  deriving  directly  from  the  p  atomic 
state  on  the  interstitial. 

The  formation  of  the  localized  levels  due  to  the  hexag- 

.al  Hg  interstitial  in  CdTe  is  shown  in  Fig.  5.  For  the 
considered  heteropolar  material  the  formation  scheme 
differs  from  that  in  a  homopolar  one  (such  as  Si),  since 
now  the  interstitial  has  both  anion  and  cation  nearest 
neighbors  (which  lowers  the  symmetry  of  the  defect  site 
to  Ch  ).  Thus,  in  panel  (a),  apart  from  s  and  p  atomic  en¬ 
ergies  on  the  Hg  interstitial,  we  now  have  both  E\  and  E2 
energies  of  the  dangling  hybrids  on  both  kinds  of  the 
INN  atoms  (i.c.,  on  Te  and  Cd  atoms,  respectively). 
When  the  orbitals  on  the  interstitial  are  allowed  to  in¬ 
teract  with  these  INN  hybrids,  in  panel  ibl  each  of  the 
E,  and  £\  levels  splits  into  a  level  of  E  symmetrv  <px,pv- 
like)  and  two  levels  of  A  ,  symmetry  I -like);  ihe  upper 
level  of  A  l  symmetry  arising  trom  the  E,  level  :s  delocal¬ 
ized  and  not  shown  in  panel  (b)  The  two  lowest  A ,  lev¬ 
els  in  panel  (b)  are  primarily  of  s  character  un  the  intersti¬ 
tial,  and  remain  so  throughout  the  entire  process  of  for¬ 
mation.  The  lower  A  {  level  arising  from  the  E:  level  is 
primarily  of  z  character  on  the  interstitial,  while  the 
upper  A  j  level  connected  by  a  dashed  line  to  the  E?  lev¬ 
el,  from  which  it  derives,  is  of  mixed  5  and  z  character  on 
the  interstitial.  The  upper  A  l  level  connected  by  a 
dashed  line  to  the  p  level  on  the  interstitial  has  a  strong  z 
character  on  the  interstitial.  \s  more  interactions  are 
added,  the  levels  shift,  but  there  are  few  additional  levels 
or  splittings.  The  additional  levels  arising  in  panel  (c), 
from  the  nclusion  of  the  covalent  interaction  between  the 
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INN  atoms,  are  hardly  localized  at  all,  and  are  not 
shown.  As  in  the  case  of  the  self-interstitial  in  Si,  when 
the  covalent  interactions  between  the  INN  and  INN 
atoms  are  added  in  panel  id),  the  strongly  localized  level 
of  z  character,  arising  from  the  p  orbitals  on  the  intersti¬ 
tial,  is  shifted  downward,  while  the  less  localized  A  j  level 
arising  from  the  £■>  level  is  shifted  upward.  Following 
the  character  of  these  strongly  mixing  state1,  the  crossing 
of  these  two  levels  is  shown  by  dashed  lines  between 
panels  tc)  and  (d).  In  panel  (d),  all  of  the  s  character  on 
the  interstitial  in  the  upper  three  A  j  states  has  shifted 
into  the  middle  one  of  these  three  states,  which  becomes 
progressively  more  delocalized  in  panels  (d)  and  (e).  As 
in  the  case  of  the  hexagonal  self-interstitial  in  Si,  the  ad¬ 
dition  of  the  promotion  energy  on  the  INN  atoms  gives 
rise  to  few  more  states.  These  includes  the  localized  /l, 
state  [at  about  —  5  cV  in  panel  te)],  which  derives  from 
the  bonding  E^  host  states  formed  due  to  the  covalent 
land  polar)  interaction  between  the  host  atoms.  The  anti¬ 
bonding  E^  host  states  give  rise  to  the  surfacelike  state 
of  A  j  symmetry  (compare  to  Sec.  Ill  C),  which  falls  be¬ 
tween  the  strongly  localized  defect  states  of  A  j  symmetry 
in  the  gap  region  in  panel  (e).  However,  now  the  charge 
on  the  defect  sue  for  this  state  is  smaller  than  for  the  Si 
self-interstitial,  ana  its  influence  on  the  localized  A 1 
states  in  the  gap  region  is  weaker;  tn  particular,  the  defect 
A  ,  state  of  strong  7  character  on  the  interstitial  remains 
above  the  defect  state  of  strong  5  character. 


IV.  CONCLUSIONS 


We  have  presented  a  simple  picture  of  the  formation 
processes  of  interstitial  defect  states,  which  allows  us  to 
interpret  the  complicated  structure  of  the  localized  levels 
in  terms  of  a  few  basic  interactions  in  the  system  consist¬ 
ing  of  the  interstitial  and  a  small  number  of  atoms  from 
its  close  vicinity.  We  have  shown  that  the  localized  states 
can  derive  either  from  atomic  states  on  the  interstitial,  or 
from  host  states  on  the  neighboring  atoms.  These  host 
states  can  be  either  simple  5/?3-hybrid  states  on  the  isolat¬ 
ed  host  atoms,  or  the  states  that  are  formed  due  to  the  co¬ 
valent  interactions  between  the  host  atoms.  The  latter 


states,  in  particular,  give  rise  to  the  mid-valence-band 
resonances  of  a  strong  hostlike  character.  We  found  that, 
while  the  structure  of  the  finally  formed  levels  is  basically 
similar  for  all  the  (respective  tetrahedral  and  hexagonal) 
interstitials  considered,  the  characters  of  some  of  the 
states  and  the  degree  of  their  localization  depend  in  an 
essential  way  on  the  ordering  (and  values)  of  the  atomic 
energies  on  the  isolated  interstitial  and  host  atoms.  For 
example,  the  Al  symmetry  tetrahedral  gap-regmn  state 
has  a  hostlike  character  for  interstitials  with  relatively 
large  values  of  s  atomic  energy  (as  for  P  in  GaP  or  the 
self-interstitial  in  Si),  while  it  may  be  strongly  defectlike 
for  interstitials  with  relatively  low  values  of  the  s  atomic 
energy  (as  for  Hg  in  CdTe/. 

We  believe  that  the  results  presented  allow  for  a  fuller 
understanding  of  the  origins  and  the  nature  of  the  local¬ 
ized  defect  states.  This  increased  understanding  may  also 
be  of  importance  in  technological  applications  involving 
defects. 
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PRESENT  STATUS  AND  FUTURE  OF  THEORETICAL  WORK 
ON  POINT  DEFECTS  AND  DIFFUSION  IN  SEMICONDUCTORS 
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A  brief  description  of  the  theoretical  methods  commonly  used 
to  investigate  point  defects  and  diffusion  is  given.  The  current 
status  of  theoretical  understanding  of  hydrogen  diffusion  in 
silicon,  self -diffusion  in  silicon,  and  Hg  and  Cd  diffusion  in 
HgCdTe  is  described,  as  examples  of  systems  where  a  basically 
coherent  picture  has  been  achieved,  where  the  many  theoretical 
contributions  have  not  succeeded  in  resolving  a  fundamental 
controversy,  and  where  theoretical  work  is  only  in  its  initial 
stages.  Several  exciting  new  directions  which  promise  expanded 
theoretical  capabilities  in  the  future  are  indicated. 
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I.  INTRODUCTION 

Veil -designed  experiments  and  theoretical  calculations  can  provide 
complementary  contributions  toward  the  ultimate  goal  of  complete  understanding  of 
point  defects  and  their  role  in  diffusion.  Although  in  some  cases,  present 
theoretical  methods  are  not  sufficiently  accurate  to  predict  which  of  several 
competing  diffusion  mechanisms  should  dominate,  they  still  offer  considerable 
guidance  in  the  search  for  a  correct  description  of  the  microscopic  atomic  motion 
during  diffusion. 

Recent  development  and  refinement  of  methods  involving  the  local  density 

approximation  (LDA)~'  to  the  density- functional  theory  for  the  exchange  and 

correlation  energy  have  made  possible  calculation  of  the  total  energies  of  the  more 

symmetric  configurations  of  a  number  of  point  defects  with  greater  quantitative 

accuracy  than  the  simpler  methods  which  had  previously  been  the  most  frequently 

used.  However,  the  amount  of  computer  time  required  by  the  new  LDA-based  methods 

can  severely  restrict  investigation  of  systems  where  the  optimal  defect 

configurations  involve  many  degrees  of  freedom  for  lattice  relaxation,  or  where 

many  different  defect  configurations  must  be  considered.  For  similar  reasons, 

LDA-based  molecular  dynamics  calculations,^  which  allow  both  defect  motion  and 

simultaneous  lattice  relaxation,  are  restricted  to  simulations  of  high  temperature 

behavior.  Fortunately,  theoretical  developments  already  in  progress  promise 

improved  accuracy  for  the  simpler  calculations  which  allow  consideration  of  a  wider 

range  of  defect  configurations,  and  may  permit  the  use  of  molecular  dynamics 

5-8 

methods  for  temperatures  which  are  currently  inaccessible,  as  well  as  improved 
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accuracy  for  the  more  elaborate  and  computationally  demanding  calculations . 

In  Section  II,  the  major  theoretical  techniques  which  have  been  used  to  study 
point  defects  and  defect-related  diffusion:  tight-binding,  Hartree-Fock-based,  and 
LDA-based  pseudopotential  and  self-consistent  Green's  function,  are  reviewed.  As 
an  example  of  successful  theoretical  characterization,  the  general  picture  of  the 
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isolated  hydrogen  interstitial  and  hydrogen  diffusion  in  silicon,  which  has  now 
been  obtained  using  a  variety  of  theoretical  approaches,  is  discussed  in  Section 
III.  Silicon  self -diffusion  is  discussed  in  Section  IV,  as  an  example  of  a  system 
where  present  Limitations  preclude  a  definitive  theoretical  determination  of  the 
point  defect  species  which  dominates  diffusion  at  high  temperature.  Section  V 
deals  with  other  diffusing  species  and  other  materials,  for  which  fewer 
calculations  have  been  done  and  many  questions  remain,  including  as  an  example  the 
current  status  of  theoretical  contributions  to  the  understanding  of  Cd  and  Hg 
diffusion  in  HgCdTe.  Section  VI  summarizes  what  we  have  learned  from  theoretical 
calculations  about  point  defects  and  diffusion,  and  indicates  some  potentially 
important  new  directions  for  the  development  of  more  convenient  and  powerful 
theoretical  tools. 

XI*  THEORETICAL  METHODS 

In  the  simplest  picture  of  thermally  activated  diffusion,  mediated  by  a  single 
point  defect  and  occurring  via  a  single  process,  the  dependence  of  the  diffusion 
rate  on  the  temperature  T  is  given  by 

diffusion  rate  -  Ce 

where  the  total  enthalpy  of  the  process  contains  a  contribution  due  to  formation  of 
the  defect,  and  a  contribution  due  to  its  migration: 

H  -  H*  .  +  H  . 

formation  migration 

Most:  of  the  theoretical  work  to  date  has  consisted  of  calculations  of  total 
energies  for  various  point  defect  configurations,  or  total  energy  surfaces  as  a 
function  of  atomic  position.  Such  work  gives  estimates  of  the  formation  and 
migration  contributions  to  the  diffusion  barrier  H,  but  it  cannot  give  absolute 
diffusion  rates,  as  it  does  not  consider  the  entropies  of  defect  formation  and 
migration,  which  enter  the  prefactor  C.  A  powerful  but  expensive  alternative 
approach,  the  molecular  dynamics  method  recently  developed  by  Car  and  Parrinello,^ 


3 


allows  simulation  of  absolute  diffusion  rates. 

The  basic  choices  which  must  be  made  in  all  of  these  calculations  are  (1)  what 
Hamiltonian  to  use  to  describe  the  electron-electron  interactions  and  the  energy  of 
the  system,  and  (2)  what  geometry  to  use  for  the  crystal,  in  order  to  best 
approximate  the  behavior  of  an  infinite  crystal  containing  a  single,  isolated 
defect . 

Empirical  tight-binding  Hamiltonians  were  originally  proposed  only  as  a  useful 

tool  for  interpolating  band  structures.  However,  the  simplicity  of  tight-binding 

calculations,  and  their  ability  to  yield  both  electronic  and  structural  energies 

simultaneously,  have  led  to  the  widespread  use  of  such  methods  to  predict 

semiquantitatively  the  electronic  and  structural  properties  of  bulk  materials, 

surfaces,  and  point  defects,  and  to  obtain  qualitative  understanding  and  physical 

insight.  Since  tight-binding  calculations  are  less  time  consuming  than  first 

principles  Hartree-Fock  or  LDA-based  pseudopotential  or  self-consistent  Green's 

function  calculations,  they  are  especially  well  suited  for  investigating  systems 

involving  many  degrees  of  freedom,  large  unit  cells,  or  complicated,  non-periodic 

17-22 

arrangements  of  atoms.  Recent  work  on  surface  reconstruction,  lattice 

23-25 

relaxation  around  tetrahedral  and  hexagonal  interstitials,  and  the  stable 

26-28 

configuration  of  small  silicon  clusters  has  shown  that  an  empirical 

tight-binding  description  of  the  total  energy  can  be  used  reliably  to  predict 

relaxation  for  systems  close  to  ideal  coordination. 

However,  in  order  to  compare  the  total  energies  of  point  defect  configurations 

with  widely  differing  atomic  coordination,  as  may  occur  when  atoms  diffuse  through 

the  crystal,  it  is  necessary  to  use  first-principles-based  expressions  for  the 

total  energy  and  interactions,  which  vary  as  a  continuous  function  of  interatomic 

distance.  Recently,  it  has  been  shown  that  various  approximations  can  transform  a 

5  7  29 

first-principles  density- functional  ’  or  muffin- tin  (MT)  calculation  into  the 
tight-binding  form,  and  a  tight-binding  approximation  to  the  LDA  pseudopotential 
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method.  valid  for  systems  with  little  or  no  charge  transfer,  has  been  successfully 

implemented  for  bulk  phases  and  small  homopolar  molecules.  Further  development 

and  testing  of  these  f irst-principles  tight-binding  methods  for  the  substantially 

ilstorted  defect  configurations,  perhaps  involving  significant  charge  transfer, 

which  may  occur  during  diffusion,  holds  the  promise  of  simpler  and  more  rapid 

calculations  of  diffusion  barriers,  and  perhaps  simulation  of  absolute  diffusion 

rates,  without  an  unacceptably  large  sacrifice  of  accuracy. 

In  contrast  to  tight-binding  methods,  which  were  originally  used  with 

empirically  fit  parameters.  Hartree-Fock  methods  were  developed  to  include 

electron-electron  exchange  as  well  as  direct  Coulomb  interactions  in  first 

principles  calculations,  with  no  adjustable  parameters.  While  accurate, 

particularly  if  expanded  and  spin-polarized  basis  sets  are  used  and  corrections  for 

electron-electron  correlation  are  included,  ab  initio  Hartree-Fock  methods  are 

computationally  demanding,  and  are  therefore  restricted  to  atomic,  molecular,  and 

small  cluster  calculations.  Use  of  larger  basis  sets  and  correlation  corrections 

increases  the  computation  time  substantially.  In  order  to  reduce  computer  costs 

and  make  these  first  principles  calculations  feasible  for  larger  systems,  certain 

small,  four-center  integrals  may  be  neglected,  resulting  in  the  partial  retention 
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of  diatomic  differential  overlap  (PRDDO)  method.  Ab  initio  Hartree-Fock  and 

PRDDO  calculations  using  minimal  basis  sets  usually  result  in  overestimates  of 

32 

potential  barriers. 

In  addition  to  the  approximate  PRDDO  method,  Hartree-Fock-based  seraiempirical 

33 

methods,  such  as  CNDO  (complete  neglect  of  differential  overlap),  MNDO  (modified 

neglect  of  diatomic  overlap),'^  and  MINDO  (modified  intermediate  neglect  of 

36 

differential  overlap) ,  have  been  developed  for  studies  involving  many 
calculations,  where  ab  initio  Hartree-Fock  methods  would  be  prohibitively 
expensive.  These  methods  use  parameters  fit  to  experimental  data  such  as  atomic 
ionization  potentials,  molecular  heats  of  formation,  and  equilibrium  bond  lengths 


in  molecules  or  crystals.  Since  single-particle  energy  differences  and 

corresponding  total  energy  differences  obtained  from  semiempirical 

Hartree-Fock-based  cluster  calculations  can  differ  substantially,  transition 

energies  cannot  be  gotten  simply  from  the  single-particle  electron  energies  in 

these  calculations,  as  they  can  in  locai-density  calculations  such  as  the 

37 

scattered-wave  Xa  (SW-Xa)"  method.  However,  when  the  parameters  are  fit  to 
structural  data  for  systems  with  similar  bonding  and  atomic  coordination,  these 
techniques  can  provide  a  good  description  of  structural  properties  and  total  energy 
changes  as  a  function  of  atomic  disp) acement. 

The  two  main  LDA-based  approaches  which  have  been  used  to  study  total  energies 
and  diffusion  are  the  self-consistent  Green's  function^  ^  and  a  priori 
pseudopotential  ^  ^  methods.  Local-density  methods  correctly  describe  bulk 

structural  properties  to  within  a  few  per  cent  of  experiment,  but  they  give  poor 
excitation  energies:  for  many  semiconductors,  the  predicted  band  gap  is  half  the 
experimental  band  gap  or  less.  For  this  reason,  the  positions  of  electronic  defect 
levels  derived  from  the  conduction  band  are  not  accurately  given,  and  total 
energies  of  defects  with  some  of  these  states  occupied  are  also  less  accurately 
given  than  bulk  total  energies.  Much  work  has  been  done  on  correcting  the  band  gap 
problem  in  bulk  calculations,^  ^  but  these  corrections  have  not  yet  been  worked 
out  for  systems  including  defects.  The  LDA  gap  has  usually  been  corrected  in 
calculations  involving  point  defects  by  the  rather  crude  procedure  of  rigidly 
shifting  the  conduction  band  and  conduction-band-derived  states  upward  to  coincide 
with  the  experimental  gap.  Estimated  error  bars  for  LDA-based  calculations  are 
about  1  eV  for  defect  formation  energies,  and  0.5  eV  for  migration  barriers.^ 

The  reliability  of  any  calculation  of  point  defect  total  energies  depends  on 
the  geometry  of  the  model  system,  as  well  as  on  the  Hamiltonian  used.  In  Green's 
function  methods,  the  distorted  defect  region  is  embedded  in  a  host  crystal 
described  by  the  perfect  crystal  Green's  function.  Although  self-consistent  LDA 
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Green's  function  calculations  have  been  carried  out  for  several  systems  (e.g. 
References  45  and  46),  these  methods  are  computationally  complex  and  very 
expensive,  particularly  if  the  defect  region  (i.  e.  the  region  which  is 

substantially  distorted  and  must  be  treated  exactly)  contains  more  than  a  few 

atoms.  In  addition,  the  role  of  the  local  environment  in  determining  the 

properties  of  the  defect  remains  obscure,  and  often  much  simpler  models  must  be 
used  to  provide  a  physical  interpretation  for  the  results  of  these  calculations. 

The  less  rigorous  but  more  physically  transparent  approach  of  the  cluster 

methods  is  to  simulate  the  infinite  crystal  by  surrounding  the  defect  with  a  finite 

number  of  atoms.  The  use  of  small  clusters  (less  than  -100  atoms)  rests  on  the 

idea  that  defect  properties  are  most  affected  by  the  local  environment  of  the 

defect.  In  cluster  calculations,  the  interactions  of  the  finite  number  of  atoms 

surrounding  the  defect  are  treated  exactly.  Rather  than  being  attached  to  the  rest 

of  the  infinite  crystal,  the  clusters  are  terminated  at  a  certain  distance  from  the 

defect,  often  by  saturating  the  dangling  bonds  at  the  surface  with  hydrogen  atoms. 

If  the  cluster  is  too  small,  interactions  of  the  defect  with  surface  states  may 

alter  the  calculated  properties.  Particularly  when  the  polarity  of  the 

host-saturator  bond  differs  from  the  polarity  of  the  bulk  host  bonds  so  that  the 

charge  transfer  at  the  surface  is  not  typical  of  the  charge  transfer  in  the  bulk, 

calculated  defect  properties  may  be  very  sensitive  to  the  host-saturator  bond 

length  used,  so  that  to  obtain  reliable  results,  care  must  be  taken  to  terminate 

47 

the  cluster  correctly. 

In  a  supercell  calculation,  the  cluster  containing  the  defect  is  not 
terminated,  but  repeated  periodically  in  space  to  form  an  infinite  crystal.  If  the 
repeated  cluster  is  too  small,  interactions  between  the  point  defects  may  produce  a 
dispersion  of  the  localized  defect  levels,  and  lead  to  calculated  total  energies 
which  are  not  the  total  energies  of  isolated  point  defects.  Therefore,  as  for 
finite  cluster  calculations,  supercell  calculations  must  be  tested  for  convergence 
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as  a  function  of  cluster  size.  In  spite  of  these  potential  pitfalls,  because  of 
their  greater  physical  transparency  and  relative  computational  simplicity,  a  priori 
pseudopotential  calculations  in  a  supercell  geometry  have  increasingly  replaced 
self 'consistent  Green's  function  calculations  as  the  most  commonly  used  LDA-based 
method. 

III.  HYDROGEN  IN  SILICON 

Not  only  is  hydrogen,  the  first  element  in  the  Periodic  Table,  a  simple  place 

to  start  in  investigating  the  properties  of  impurities  in  semiconductors,  it  is  an 

impurity  with  a  major  technological  significance.  Hydrogen  can  passivate  both 

donor  and  acceptor  dopants,  as  well  as  lattice  defects  (see,  e.  g. ,  References 

48-51,  and  Reference  52  for  a  review).  When  added  in  sufficient  concentrations, 

hydrogen  can  form  a  molecular  species,  and  inhibit‘d  or  contribute  to"^  the 

formation  of  extended  defects.  Since  hydrogen  interacts  strongly  with  so  many  of 

the  impurities  and  lattice  imperfections  that  may  be  found  in  crystals,  the 

fundamental  diffusion  behavior  of  hydrogen  in  a  perfect  semiconductor  crystal  has 

been  difficult  to  determine  experimentally.  Both  different  concentrations  of 

trapping  sites  for  hydrogen,  and  the  space  charge  fields  due  to  variations  in 

dopant  or  hydrogen  concentration,  which  may  affect  the  motion  of  charged  hydrogen, 

have  resulted  in  the  wide  variety  of  results  reported  for  hydrogen  diffusion  at 
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temperatures  below  800  K. 

In  order  to  avoid  the  complex  details  of  the  interactions  of  hydrogen  with 
various  dopants,  where  much  still  remains  controversial,  we  will  discuss  here  only 
the  behavior  of  an  isolated  hydrogen  atom  in  a  perfect  silicon  crystal.  The  same 
basic  picture  of  the  total  energy  for  this  system  as  a  function  of  hydrogen 

position  is  now  supported  by  calculations  using  a  variety  of  theoretical 

techniques. 

In  order  to  describe  how  hydrogen  moves  through  the  crystal,  the  stable  and 
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metastable  positions  where  hydrogen  prefers  to  sit  must  first  be  identified.  In 

contrast  to  kinetic  behavior,  the  location  of  these  sites  is  independent  of  the 

mass,  and  depends  only  on  the  chemical  nature  and  bonding  properties  of  hydrogen. 
Therefore,  the  stable  and  metastable  positions  in  the  lattice  should  be  the  same 
for  hydrogen  as  for  muonium,  which  may  be  considered  to  be  a  light  isotope  of 
hydrogen.  There  is  a  substantial  amount  of  experimental  data  on  the 
characteristics  of  muonium  in  semiconductors,  which  is  reviewed  in  Reference  59. 
In  muon  spin  rotation  experiments,  a  muonium  state  with  an  isotropic  hyperfine 

interaction  (’'normal"  muonium)  and  a  muonium  state  with  a  hyperfine  interaction 
that  has  <111>  axial  symmetry  ("anomalous"  muonium)  are  observed  in  silicon, 
diamond,  and  several  other  semiconductors  with  diamond  or  zincblende  lattice 


structure . 

Early  theoretical  work,  including  extended  Huckel  calculations,^  Green's 

function  calculations,^  MNDO  calculations,^  and  CNDO  calculations^  for  hydrogen 

64 

in  silicon,  Hartree-Fock  calculations  tor  muonium  in  diamond,  and  later 
Hartree-Fock  and  PRDDO  calculations  for  muonium  and  hydrogen  in  diamond*^  and  in 
silicon, ^  considered  both  high-electron-density  sites,  on  or  near  the  chains  of 
bonded  host  atoms,  and  low-electron-density  sites,  in  the  open  spaces  between  the 
bonding  chains.  These  calculations  usually  included  little  or  no  relaxation  of  the 
host  lattice,  and  tended  to  find  low-electron-density  sites  as  the  energetically 
favored  positions  for  hydrogen  or  muonium  in  the  lattice.  (Figure  1  shows  the 
position  of  the  low-electron-density  tetrahedral  and  hexagonal  sites  and  the 
high-electron-density  bond-centered  site  in  the  lattice.  For  clarity,  not  all  the 


positions  indicated  by  various  papcrr 


as  actua^ 


m  the  potential  energy 


surface  are  shown  in  this  figure.)  Green's  function  calculations  for  a  muon 
assumed  to  be  at  the  tetrahedral  interstitial  site  in  silicon^  obtained  a  value 


very  close  to  the  experimental  value  for  the  hyperfine  coupling  constant  of  normal 
muonium.  While  the  earliest,  simple  treatments  were  later  superseded  by  more  exact 


E-iO 
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calculations,  and  some  of  the  early  calculations  were  later  refined,  as 

understanding  of  how  to  eliminate  surface-defect  interactions  due  to  cluster  size 

and  surface  termination  progressed,  the  conclusion  that  low-electron-density  sites 

at  or  near  the  tetrahedral  position  were  favored  energetically  for  mi  ' ^al  lattice 

relaxation  remained.  (See,  for  example,  the  later  Hartree-Fock  and  FRDDO 

65 

calculations  for  hydrogen  and  muonium  in  diamond,  which  corrected  the  earlier 
treatment  of  bond  saturation  at  the  surface,  ^  shown  in  Reference  68  to  lead  to 
cluster  size  dependent  results.)  More  recent  calculations^^  ^  do  not  contradict 
the  identification  of  normal  muonium  in  silicon  with  a  muonium  traveling  (or 
tunneling)  from  one  site  to  another  in  the  low-electron-density  channels,  where  the 
total  energy  surface  is  relatively  flat,  but  trapped  in  these  channels  by  the  lack 
of  time  for  the  large  lattice  relaxations  required  to  stabilize  the  muonium  at 
other  positions. ^  A  similar  motion  of  hydrogen  through  the  low-electron-density 
channels  in  silicon  has  been  observed  in  high  temperature  molecular  dynamics 
calculations . ^ 

The  bond-centered  position  was  first  suggested  as  the  site  for  anomalous 

72 

muonium  by  Cox  and  Symons,  although  Fisch  and  Licciardello  had  previously 

proposed  this  site  as  a  likely  position  for  hydrogen  in  amorphous  silicon,  where  a 

range  of  silicon-silicon  bond  lengths,  including  some  of  the  proper  length  to 

73 

accommodate  an  added  hydrogen  atom,  is  expected.  According  to  the 

three-center-bond  model  described  in  Reference  73,  the  symmetric  combination  of  the 
sp  hybrids  on  the  two  silicon  neighbors  which  point  toward  the  hydrogen  interacts 
with  the  hydrogen  s  orbital  to  form  a  totally  symmetric,  bonding  wave  function  and 
an  antibonding  wave  function,  while  the  antisymmetric  combination  of  the  silicon 
hybrids  does  not  interact  with  the  hydrogen  s  orbital,  and  remains  a  non-bonding 
wave  function  for  the  three-center  geometry.  Since  each  silicon  contributes  one 
electron  to  fill  these  electronic  states,  for  positively  ionized  hydrogen  H  ,  which 
has  no  electron  to  contribute,  only  th  three-center  bonding  state  should  be 


10 


filled.  The  chree-center  geometry  should  therefore  be  stable  for  H  ,  if  the 

nearest  neighbor  silicons  are  allowed  to  relax  to  their  preferred  distance  from  the 

hydrogen.  For  neutral  hydrogen,  one  electron  must  go  in  the  non-bonding  state,  so 

the  three-center-bond  configuration  is  less  likely  to  be  stable,  and  for  negatively 

charged  H  ,  with  the  non-bonding  state  completely  filled,  this  configuration  is 

even  less  favorable.  The  identification  of  the  bond-centered  site,  stabilized  by  a 

large  relaxation  of  the  neighboring  silicon  atoms,  as  the  preferred  position  ot 

muonium  or  hydrogen  in  silicon  was  confirmed  initially  by  the  PRDDO 

calculations,^  ^  and  subsequently  by  MINDQ,^  ^  MNDO,^  and  LDA 

pseudopotential  *  calculations.  The  general  picture  emerging  from  these 

calculations  is  that  neutral  hydrogen  and  H+  prefer  to  sit  at  the  bond-centered 

site,  while  H  prefers  a  low-electron-density  site.  The  theoretical  identification 

74-75 

of  the  bond-centered  site  as  the  position  of  anomalous  muonium  was  also 

verified  subsequently  by  a  muon  spin  resonance  experiment  showing  the  presence  of 
two  equivalent  silicon  nearest  neighbors,  each  bonded  to  three  other  silicon 
atoms,  and  an  experimental  observation  of  a  hydrogen  center  in  silicon  which  has 
a  hyperfine  structure  similar  to  that  of  anomalous  muonium,  and  should  therefore 


also  be  located  at  the  bond-centered  site,  was  reported. 
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The  total  energy  surface  as  a  function  of  hydrogen  position  obtained  from  LDA 

79 - 80  84 

pseudopotential  calculations  ’  has  been  used  to  calculate  diffusion 

coefficients  for  hydrogen  in  silicon,  assuming  the  silicon  lattice  is  allowed  to 

85 

relax  completely  for  each  hydrogen  position.  In  this  work,  it  is  shown  that  the 

most  energetically  favorable  path  for  H+  (assuming  complete  lattice  relaxation) 

lets  the  hydrogen  wind  its  way  down  the  bonding  chains  between  bond- centered  sites 

without  ever  going  out  into  the  low-electron-density  channels,  while  the  low  and 

84 

high-electron  density  paths  are  similar  in  energy  for  neutral  hydrogen.  However, 
it  is  unlikely  that  the  much  more  massive  silicon  atoms  have  time  to  relax 
substantially  in  response  to  the  motion  of  the  hydrogen,  as  was  pointed  out,  for 
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example,  in  Reference  86,  so  that  the  motion  of  hydrogen  between  the  bond-centered 

sites  is  likely  to  be  determined  primarily  by  the  ability  of  tne  silicon  lattice  to 

assume  a  position  favorable  for  the  new  bond-centered  location.  Molecular  dynamics 

simulations^  of  the  high  temperature  diffusion  of  H+  show  that  the  hydrogen  most 

frequently  jumps  between  bond-centered  sites  by  passing  through  the 

low-electron-density  hexagonal  and/or  tetrahedral  sites,  but  it  can  also  follow  a 

path  lying  completely  in  the  low-electron- density  region.  Due  to  the  large  energy 

barriers  which  oppose  insertion  of  the  hydrogen  into  a  new  bond-centered  site 

before  the  neighboring  silicon  atoms  have  had  time  to  relax,  tunneling  effects 

(which  have  not  been  included  in  any  of  the  work  on  hydrogen  diffusion  through  a 

perfect  crystal)  are  not  expected  to  change  the  picture  substantially,  although 

tunneling  of  hydrogen  around  a  beryllium  atom  to  which  it  is  bound  in  a  silicon 

87 

crystal  has  been  both  observed  and  explained  theoretically,  since  in  this  case  an 

energetically  favorable  path  for  the  motion  which  does  not  involve  substantial 

88 

lattice  relaxation  has  been  found. 

In  summary,  although  details  of  the  potential  energy  surface  calculated  by 
different  theoretical  methods  may  differ  and  some  questions  remain,  the  general 
picture  for  hydrogen  in  a  perfect  silicon  crystal  is  remarkably  clear.  One  major 
lesson  that  has  emerged  from  this  work,  which  was  not  initially  obvious,  is  that 
thorough  consideration  of  lattice  relaxation  is  essential,  and  highly  distorted 
configurations  of  point  defects  can  turn  out  to  be  the  most  energetically 
favorable.  Unfortunately,  no  similarly  clear  conclusions  about  which  of  the 
relevant  processes  is  dominant  at  high  temperature  can  as  yet  be  drawn  from  the 
theoretical  work  on  self -diffusion  in  silicon. 

IV.  SILICON  SELF-DIFFUSION 

The  variety  of  processes  and  defect  configurations  which  may  be  important  in 
self -diffusion  is  larger  than  for  diffusion  of  hydrogen  in  a  perfect  crystal. 
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Diffusion  may  occur  by  motion  of  vacancies  as  well  as  interstitials,  or  even  by 

direct  exchange  of  atoms  on  the  host  lattice,  in  the  absence  of  point  defects.  A 

diffusing  self  -  interstitial  may  try  to  share  a  lattice  site  with  a  host  atom,  in 

any  of  several  split- interstitial  configurations.  A  large  amount  of  experimental 

work  has  been  done  on  diffusion  in  silicon:  for  surveys  of  the  more  recent  work  and 

references  to  earlier  reviews,  see  References  89  and  90. 

Vacancies  were  the  first  native  point  defects  in  silicon  to  be  extensively 

39-40  91  92-93 

studied  both  theoretically  *  and  experimentally.  A  coherent  and 

reasonably  complete  picture  of  the  isolated  vacancy  has  resulted.  As  predicted 

91  94-95  96 

theoretically,  ’  the  vacancy  was  later  experimentally  shown  to  be  an 

97  + 

Anderson  negative-U  system:  in  other  words,  the  singly  charged  V  state  was  shown 

0  -H- 

to  be  metastable,  decaying  either  into  V  or  V  depending  on  the  Fermi  energy. 

98-101 

The  neutral  vacancy  formation  enthalpy  of  about  3.5  eV,  and  migration 

102  0  {  { 

enthalpies  for  the  direct  path  along  a  bond  of  0.3  eV  for  V  and  0.4  eV  for  V  , 

calculated  by  LDA  methods,  and  allowing  the  lattice  to  relax  completely  for  each 

position  of  the  traveling  atom,  agree  well  with  the  experimental  values  of 
103-104  92-°3 

3.4eV,  0.33  eV,  and  0.45  eV,  "  respectively.  Further  LDA  calculations 

show  that  migration  by  a  distorted  path,  preserving  as  many  bonds  to  the  traveling 
atom  as  possible,  seems  to  be  more  energetically  favorable  than  vacancy  migration 
by  the  direct  path,  along  the  bond.^^ 

Although  self -diffusion  in  silicon  appears  to  be  mediated  by  vacancies  at  low 
89 

temperature,  the  identification  of  the  dominant  high  temperature  diffusion 
mechanism  remains  controversial.  For  an  interpretation  of  the  experimental 
evidence  which  supports  the  conclusion  that  vacancies  and  not  interstitials 
dominate  diffusion  at  high  temperature,  see  Reference  106.  LDA  simulations  of  high 
temperature  vacancy  diffusion  show  strong  dynamic  effects,  which  may  help  explain 
the  large  pre- exponential  factors  seen  in  high  temperature  diffusion. 

In  contrast  with  vacancies,  isolated  self -  interstitials  have  never  been 


S-14 


observed  in  silicon.  Even  after  radiation  at  4.2  K,  the  displaced 

self- interstitial  appears  to  migrate  rapidly  until  it  displaces  a  substitutional 

impurity  (or  is  otherwise  trapped) ,  resulting  in  the  observation  of  impurity 

interstitials  or  interstitial-related  defect  complexes,  but  no  isolated 

108 

self  -  interstitials .  Early  extended  Hiickel  calculations  for  the  carbon 

interstitial  in  diamond, which  were  used  to  suggest  similar  behavior  for  the 

112 

silicon  self - interstitial ,  suggested  that  athermal  migration  might  be  possible, 

since  the  predicted  minimum- energy  configuration  changed  for  different  interstitial 
charge  states.  These  calculations  also  indicated  an  apparent  preference  of  the 
reutral  interstitial  for  a  bond-centered  or  ^plit  configuration  rather  than  a 
low-electron-density  tetrahedral  or  hexagonal  interstitial  site. 

More  accurate  and  elaborate  LDA-based  calculation  were  later  performed  for  the 


tetrahedral  and  hexagonal  configurations  of  the  silicon  self -  interstitial , 
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as  well  as  for  several  varieties  of  the  split  and  bond-centered 
98  100,116  117  m  •  M  Pi  a  v>n  1  t.tU  i  mU  a.tAVn  VO  i  PPi  /^1  I  1  f*  h  A  t-  >-  /-a 


configurations , 


shown  in  Figure  1,  which  were  more  difficult  to  treat 


accurately,  due  to  their  lower  symmetry  and  larger  lattice  relaxations.  In  all 

i  | 

these  calculations,  the  doubly  ionized  self- interstitial  Si  was  found  to  prefer 
the  tetrahedral  site,  where  the  electron  density  is  the  lowest,  while  the  neutral 
interstitial  was  found  to  prefer  sites  of  higher  electron  density,  such  as  the 
hexagonal,  bond- centered,  and  split  interstitial  sites.  Of  all  the  configurations 
considered  in  these  static  calculations,  the  <110>  split  configuration,  shown  in 
Figure  1  as  the  X  configuration,  was  found  to  have  the  lowest  energy  for  the 
neutral  interstitial . However,  the  self -interstitial  was  found  to  have  similar 
total  energies  in  several  configurations.  This  should  result  in  similar  activation 
enthalpies  for  several  migration  paths  and  large  activation  entropies.  The  results 
also  suggest  that  there  may  be  many  paths  available  for  athermal  diffusion  due  to 
charge  carrier  capture.  LDA-based  molecular  dynamics  calculations  in  progress 


have  so  far  confirmed  that  the  neutral  X  interstitial  is  stable  and  that  there  are 
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large  flat  regions  in  the  total  energy  surface  as  a  function  of  interstitial 

position,  which  should  give  rise  to  high  entropies  of  formation  for  silicon 

118 

self- interstitials . 

In  summary,  both  vacancies  and  interstitials  have  been  found  to  have  formation 

enthalpies  of  roughly  4  eV  and  migration  enthalpies  of  roughly  0.4  eV,  giving  rise 

89 

to  activation  enthalpies  for  diffusion  in  the  experimental  range  of  4-5  eV.  The 

large  activation  enthalpies  observed  at  high  temperature,  where  defects  must  be 

created  thermally,  and  the  rapid  diffusion  seen  in  low  temperature  radiation 

experiments,  can  therefore  be  explained  by  the  same  point  defect  mechanisms. 

However,  since  the  calculated  activation  enthalpies  for  vacancies  and  interstitials 

are  similar,  it  is  not  possible  to  resolve  the  controversy  over  the  dominant 

mechanism  for  silicon  self -diffusion  at  high  temperature  on  the  basis  of  presently 

available  theoretical  results.  To  make  things  still  more  complicated,  the  direct 

exchange  mechanism  proposed  by  Pandey,  which  takes  place  without  point  defects, 

119 

also  has  an  activation  enthalpy  in  the  experimental  range. 

The  situation  may  become  somewhat  clearer  in  the  near  future,  with  the 

publication  of  the  results  of  molecular  dynamics  simulations  of  high  temperature 

diffusion  mediated  by  vacancies'^  and  interstitials .  In  the  meantime,  static 

LDA-based  calculations  of  the  effects  of  pressure  on  the  activation  energies  for 

diffusion  by  vacancy,  interstitial,  and  concerted  exchange  mechanisms  have 

suggested  that  an  experimental  investigation  of  diffusion  under  different  pressures 

120 

be  undertaken  to  determine  the  dominant  mechanism. 


V.  OTHER  SYSTEMS:  Hg,  Cd  Te 

x-X  x 

Since  the  same  native  defects  which  dominate  self -diffusion  may  also  play  a 
role  in  the  diffusion  of  dopants,  work  on  dopant  diffusion  in  silicon  has  been 
closely  tied  to  work  on  self -diffusion.  Because  of  its  intrinsic  technological 


importance  as  well  as  its  relevance  for  self-diffusion,  there  has  been  much 
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experimental  work  on  dopant  diffusion  in  silicon  (surveyed  in  References  89  and  90, 

along  with  some  discussion  of  the  points  which  have  remained  controversial  for 

dopant  diffusion  as  well  as  self -diffusion) ,  as  well  as  a  number  of  theoretical 

calculations.  Reference  115,  for  example,  shows  that  carrier-capture -enhanced 

diffusion  is  possible  for  both  aluminum  and  self-interstitials  in  silicon.  Other 

LDA-based  calculations^^  ’  indicate  that  diffusion  for  boron,  phosphorus, 

arsenic,  and  antimony  in  silicon  may  proceed  down  interstitial  channels  or  via  a 

119 

vacancy-assisted  process,  while  the  concerted  exchange  mechanism,  which  involves 

no  point  defects,  appears  to  be  less  important.  Initial  LDA-based  calculations  of 

the  total  energy  of  chalcogen  impurities  at  various  sites  in  silicon,  and 

122-123 

preliminary  applications  for  diffusion,  have  been  performed.  Investigation 

of  the  rich  variety  of  defect  reactions  which  may  be  involved  in  the  diffusion  of 

oxygen,  carbon,  and  other  impurities  in  silicon  by  empirical  classical 

potentials ,  MNDO ,  MINDO,^^  ^  and  other  methods  has  only  begun. 

Even  less  is  known  about  impurity  and  self -diffusion  in  most  other  systems. 

In  addition  to  the  early  semiempirical  work  on  diffusion  in  diamond,  including  the 

extended  Huckel  calculations  mentioned  above,  which  suggested  that  the  carbon 

self - interstitial  prefers  a  bond-centered  or  split  configuration  rather  than  a 

tetrahedral  or  hexagonal  configuration, considerably  more  accurate  recent 

LDA-based  calculations  have  indicated  that  vacancies  rather  than  interstitials  or 

direct  exchange  dominate  self -diffusion  in  diamond  for  all  positions  of  the  Fermi 
128  -  9 

level.  Similar  calculations  of  defect  formation  energies  in  SiC  show  that  the 

antisite  defect  is  favored  in  C-rich  material,  while  for  Si-rich  material,  the 

antisite  defect  Sir  is  favored  in  n-type  crystals  and  the  carbon  vacancy  is  favored 

129-130 

in  p-type  crystals.  Recent  calculations  of  the  total  energies  of  various 

131-132 

defects  in  GaAs  form  an  initial  starting  point  for  an  LDA-based 

133 

investigation  of  diffusion  in  this  material.  Other  LDA-based  calculations  show 
that  a  vacancy  mechanism  appears  to  be  favored  for  self -diffusion  in  graphite, 
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although  they  cannot  rule  out  the  possibility  of  a  more  complicated  direct  exchange 
process  than  the  one  considered  (as  for  ail  static  calculations,  future 
consideration  of  other,  more  complex  diffusion  processes,  perhaps  involving  more 
highly  distorted  point  defect  configurations,  may  show  these  to  be  favored  over  the 
processes  previously  considered) . 

Behavior  observed  in  binary  and  ternary  compounds  involves  defects  which  may 
contain  anion  or  cation  atoms,  and  may  be  situated  at  a  larger  number  of 
distinguishable  sites  than  are  found  in  the  diamond  lattice  of  a  homopolar  crystal. 
For  example,  two  types  of  tetrahedral  interstitial  site  occur  in  the  zincblende 
crystal:  sites  with  cation  and  with  anion  nearest  neighbors,  which  we  may  here  call 
T1  and  T2  sites,  respectively.  The  number  of  available  defects  and  defect 
complexes  in  the  these  systems  and  the  variety  of  configurations  which  may  be 
important  for  diffusion  is  so  great  that  initial  theoretical  calculations  can  only 
hope  to  give  some  additional  direction  for  further  experimental  and  theoretical 
work. 

The  current  state  of  understanding  of  Cd  and  Hg  diffusion  in  CdTe  and 

Hg-  Cd ^Te  alloys  may  be  considered  to  be  illustrative  of  the  situation  for  many  of 

the  III-V  and  II-VI  semiconductors,  where  a  reproducible  experimental 

characterization  of  the  diffusion  does  not  always  bring  with  it  a  coherent 

understanding  of  the  microscopic  process,  or  even  which  point  defects  are 

responsible  for  the  diffusion,  and  how  they  contribute  to  observed  electrical 

behavior  and  other  properties  of  the  material.  Recent  tracer  diffusion 
133  1 AO 

results  on  self -diffusion  of  Cd,  Hg,  and  Te  in  Hg^  gCd^  ^Te  have  produced  a 

reasonably  reproducible  experimental  characterization,  despite  considerable 
inconsistency  in  the  early  experimental  work  (reviewed  in  References  134  and  135). 
Mercury  profiles,  which  are  unchanged  by  preannealing,  may  be  fitted  by  a  sum  of 
two  complementary  error  function  (ERFC)  curves.  The  slower  component  appears  to  be 
independent  of  the  mercury  partial  pressure,  p  ,  while  the  diffusion  coefficient 
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of  the  fast  component  is  proportional  to  p^  at  high  mercury  pressures  and  to  p^ 
at  lou  „rcury  pressures . (Brown  and  Willoughby  ’  see  additional, 
faster  components  after  diffusion  anneals  at  225°C,  instead  of  the  slow  component 
seen  in  References  135-138  and  140,  which  may  have  been  too  slow  to  observe  at 
their  diffusion  temperature  and  annealing  time.)  Two  components  are  also  seen  for 
Cd  diffusion  (except  for  Reference  139,  which  reports  only  one  component).  At 
350°C  and  above,  Cd  diffusion  appears  to  be  independent  of  the  mercury  partial 

pressure  (which  determines  the  partial  pressures  of  the  other  components  as 

, . 142.  ^  .  135-139 

well  ),  except  for  an  increase  near  the  maximum  value  of  p^. 

The  behavior  of  Cd  in  Hgn  0Cdn  0Te  is  similar  to  Cd  self -diffusion  behavior 
observed  in  CdTe  (reviewed  in  Reference  144).  Aside  from  the  small  increase  in  Cd 
diffusion  observed  for  CdTe  under  higher  Cd  pressures  (corresponding  to  the  similar 
increase  observed  in  HgQ  gCd^  ^Te)  ,  which  has  been  attributed  to  the  increasing 
importance  of  Cd  interstitials,  the  dominant  mechanisms  suggested  for  Cd  diffusion 
in  CdTe  include  a  pressure- independent  exchange  process  which  does  not  involve 
point  defects,  or  a  process  involving  a  neutral  defect  complex,  perhaps  consisting 
of  a  Cd  vacancy  and  a  Cd  interstitial .  Shaw^^  suggests  that  complexes  of 

native  defects  could  account  for  Cd  diffusion  in  Hg^  gCdQ  2Te*  TanS  anci 
Stevenson  suggest  that  a  complex  mechanism  involving  both  a  Cd  interstitial  and 
a  cation  vacancy  is  responsible  for  the  pressure  independence  of  Cd  diffusion, 
while  vacancy  and  interstitial  mechanisms  operating  in  parallel  are  responsible  for 
the  pressure  dependence  observed  for  the  fast  component  of  Hg  diffusion  in 
Hg()  8Cc*0  2^6,  However,  the  Tang  and  Stevenson  model  of  vacancy- interstitial 
diffusion  in  series  for  Cd  diffusion  does  not  explain  why  the  interstitial  step 
must  be  followed  immediately  by  the  vacancy  step,  instead  of  occurring  repeatedly 


as  a  purely  interstitial  kick-out  diffusion  process,  in  parallel  with  any  existing 

135 

vacancy-mediated  processes.  Tang  and  Stevenson's  belief  that  the  fast  component 


seen  for  mercury  represents  true  serf-diffusion  in  the  bulk,  rather  than  some  sort 
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of  short-circuit  or  pipe  diffusion,  has  also  been  questioned. 

The  earliest  assistance  for  experimental  interpretation  and  planning  in  HgCdTe 
was  obtained  from  semiempirical  estimates  of  bond  energies  and  point  defect 
formation  energies.  These  estimates  have  now  been  followed  by  tight-binding 
^  ^  -  * - 1  - -  -  -  - 1  - ' — 1 - .*^-25,147  ancj  interstitials  in 


calculations  of  the  total  energies  of  isolated  vacancies 
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various  configurations,  ’  and  subsequently  linear  muffin- tin  orbital 

148 

calculations  using  the  atomic  spheres  approximation  (LMTO-ASA)  for  interstitials 

in  the  two  tetrahedral  positions  and  antisite  defects.  Lattice  relaxation  was 

included  in  some  of  the  interstitial  calculations  (References  23-25);  however,  a 

complete  investigation  including  highly  distorted  interstitial  configurations  has 

not  yet  been  done,  and  no  defect  complexes  have  yet  been  considered. 

Specific  conclusions  which  may  be  drawn  from  this  work  include  the  suggestion 

that  Te  antisite  defects  may  be  more  important  than  previously  assumed:  these 

defects  appear  to  have  a  favorable  total  energy,  and  the  creation  of  these 

defects  may  be  promoted  in  Te-rich  material  by  a  strain-mediated  interaction 

25 

between  Te  interstitials  and  cation  vacancies.  The  attraction  of  interstitials, 

which  have  an  outward  nearest-neighbor  relaxation,  to  vacancies,  which  have  an 

inward  nearest-neighbor  relaxation,  should  be  particularly  strong  for  Te,  since  the 

calculated  relaxation  for  the  preferred  tetrahedral  interstitial  position  was 

25 

larger  for  Te  than  for  all  the  other  native  or  dopant  interstitials  considered. 

The  large  strain  energy  associated  with  Te  interstitials  may  also  help  to  explain 

25 

the  tendency  for  excess  Te  to  condense  into  second  phase  inclusions. 

Interstitial  diffusion  down  the  open  channels  between  tetrahedral  sites, 

n  /.  6 

passing  through  the  hexagonal  sites,  was  initially  considered  by  comparing  total 

energies  at  the  two  tetrahedral  sites,  as  the  original  tight-binding  treatment  did 

not  permit  a  reliable  direct  comparison  of  total  energies  at  interstitial  sites 

6*8  26-28 

with  different  coordination.  ’  Estimated  diffusion  barriers  obtained  in 

this  way  indicated  that  direct  interstitial  diffusion  down  the  open  channels 
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appeared  to  be  energetically  favorable  for  Hgt  but  unfavorable  for  Cd,  which  must 

therefore  diffuse  by  some  more  complicated  mechanism. Subsequent  calculations 

showing  that  large  strains  (and  strain  energies)  are  associated  with  Cd 

interstitials  at  T1  and  hexagonal  sites^  may  help  explain  the  apparent 

unfavorability^^  ’  of  Cd  interstitial  diffusion  down  the  open  channels  between 

preferred  tetrahedral  (T2)  sites.  For  Hg  diffusion,  where  a  simple  interstitial 

mechanism  appears  to  be  energetically  favorable ,  the  possibility  that 

146 

diffusion  could  be  enhanced  by  charge  carrier  capture  was  suggested,  since  the 
preferred  charge  state  changes  with  position, going  from  neutral  in  the  T1 
position  to  doubly  charged  Hg  in  the  preferred  tetrahedral  (T2)  and  hexagonal 
positions . 

Obviously  much  more  work,  both  theoretical  and  experimental,  will  need  to  be 

done  before  the  atomic  processes  which  occur  during  Hg  and  Cd  diffusion  in  CdTe  and 

Hg^  ^Cd^Te  are  fully  understood.  All  this  work  will  be  complicated  by  the 

interactions  of  native  point  defects  with  impurities,  and  with  any  strains  which 

may  be  present  in  the  lattice  as  a  result  of  nearby  extended  defects,  surfaces,  or 

interfaces.  For  example,  a  strong  dependence  of  point  defect  formation  energies  on 

these  strains  in  the  lattice  could  result  in  strongly  inhomogeneous  point  defect 

concentrations  and  diffusion  constants,  perhaps  producing  effects  such  as  the 

151-152 

depth-dependent  diffusion  recently  observed  in  HgCdTe  superlattices. 


VI.  CONCLUSIONS 

The  characteristics,  advantages,  and  shortcomings  of  the  basic  theoretical 
methods  used  to  investigate  point  defects  and  diffusion  in  semiconductors  have  been 
briefly  reviewed.  While  these  methods  differ  widely  in  accuracy  and  computational 
expense,  all  use  approximations  which  cause  errors  in  point  defect  total  energies 
and  diffusion  barriers  that  cannot  be  precisely  estimated,  except  by  appeal  to  test 
calculations  on  similar  systems.  Both  the  development  of  tight-binding 


20 
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calculations  based  on  first-principles  Hamiltonians,  for  the  initial  exploration  of 
systems  where  many  defects  or  many  configurations  must  be  considered,  and 
corrections  leading  to  an  improved  band  gap  for  the  more  accurate  and  expensive 
LDA-based  defect  calculations,  should  lead  to  expanded  theoretical  capabilities  in 
the  foreseeable  future. 

The  history  and  current  status  of  theoretical  work  on  hydrogen  diffusion  in 
silicon,  silicon  self -diffusion,  and  diffusion  of  Hg  and  Cd  in  HgCdTe  have  also 
been  described,  as  examples  of  a  system  where  a  basically  coherent  picture  has  been 
achieved,  a  system  where  the  many  theoretical  contributions  have  still  not  resolved 
a  fundamental  controversy,  and  a  system  where  only  the  very  first  steps  have  been 
taken  toward  a  theoretical  understanding  of  diffusion.  Aside  from  the  specifics  of 
what  has  been  learned  for  each  system,  the  theoretical  calculations  on  hydrogen  and 
self - interstitials  in  silicon  have  contributed  to  the  growing  realization  of  the 
importance  of  point  defect  configurations  involving  considerable  relaxation  of 
neighboring  atoms.  The  realization  that  the  preferred  configuration  of  a  defect  is 
not  always  the  most  obvious  one  has  made  the  value  of  the  new  molecular  dynamics 
methods,  which  allow  the  defect  and  the  surrounding  lattice  to  relax  freely  and 
simultaneously,  very  apparent.  This,  together  with  the  new  understanding  of 
dynamic  effects,  unexplained  by  static  total  energy  calculations,  which  is  being 
reached  by  molecular  dynamics  calculations  for  hydrogen  and  self -diffusion  in 
silicon,  makes  the  advent  of  these  methods  the  third  exciting  new  development,  from 
which  much  can  be  expected  in  the  future. 
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FIGURE  CAPTION 


Fig,_  1 .  The  tetrahedral  (I),  hexagonal  (H) ,  and  bond-centered  (B)  interstitial 
positions  in  the  diamond  or  zincblende  lattice.  Two  types  of  split  interstitial 
are  shown,  with  the  two  atoms  which  share  one  atomic  site  in  the  host  lattice 
relaxed  along  the  <110>  (X)  and  <001>  (S)  directions. 


